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Abstract- Direct imaging and spectroscopy is the likely means by which we will someday 
identify, confirm, and characterize an Earth-like planet around a nearby Sun-like star. This Chapter 
summarizes the current state of knowledge regarding discovering and characterizing exoplanets 
by direct imaging and spectroscopy. We detail instruments and software needed for direct imaging 
detections and summarize the current inventory of confirmed and candidate directly-imaged 
exoplanets. Direct imaging and spectroscopy in the past decade has provided key insights into 
jovian planet atmospheres, probed the demographics of the outskirts of planetary systems, and 
shed light on gas giant planet formation. We forecast the new tools and future facilities on the 
ground and in space that will enhance our capabilities for exoplanet imaging and will likely image 
habitable zone rocky planets around the nearest stars. 


1. Introduction 


Identifying, confirming, and characterizing an Earth-like 
planet around a nearby Sun-like star is a key goal of exo- 
planet science. In the coming decades, direct imaging and 
spectroscopy are the likely means by which this goal will 
be achieved. 

Over 5000 extrasolar planets and candidates have been 
detected through indirect means. In contrast, 20—25 exo- 


planets have been directly imaged (e.g. |Marois et al.|2008b 
Lagrange et al.|2010| see Fig.[1]. Despite this relatively low 


yield of discoveries, direct imaging comprises a large frac- 
tion of the known exoplanet population amenable to atmo- 
spheric characterization, since the method provides photons 
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from the planets themselves instead of just inferring their 
presence (e.g. Currie et al.|2011] Konopacky et al.[2013). 

Almost all of the directly imaged planets known thus 
far are young, self-luminous gas giants detected in thermal 
emission. They represent the extremes of planet formation 
Vigan et al.|2021): planets with masses 
of 2 Mj or more orbiting at moderate to wide separations 
(~ 10-250 au). A subset of directly detected planets appear 
to be in or just finished with the final stages of assembly 
(“protoplanets”) (Keppler et al.|2018). The direct imaging 
field currently also provides key insights into jovian planet 
formation and the demographics of the outskirts of plane- 
tary systems. 

Remarkably, the direct imaging method enables the ac- 
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Fig. 1.— Demographics of planetary-mass companions detected 
with various methods. Bold symbols denote the 175 compan- 
ions with spectroscopic measurements constraining planet atmo- 
spheres: directly imaged companions comprise a substantial frac- 
tion of this population. Figure courtesy of Dmitry Savransky, us- 
ing data from the NASA Exoplanet Archive. 


quisition of hundreds to thousands of spectroscopic data- 
points on exoplanets in a few hours of telescope time. This 
wealth of information can critically constrain individual 
planet atmospheric properties — e.g. temperature, clouds, 
chemistry, and gravity — as well as the atmospheric evolu- 
tion of gas giants as a population. 

The first direct images and spectra of planet candidates 
relied on facility, general-use adaptive optics (AO) systems 
to deblur starlight or space telescopes (i.e. the Hubble 
Space Telescope) coupled with simple coronagraphs (e.g. 
Kalas et al.[2008). The past decade has 
seen the development, demonstration, and honing of dedi- 
cated extreme AO systems coupled with advanced coron- 
agraphs and sophisticated post-processing methods, allow- 
ing the detection of planets that are fainter and less mas- 
sive and/or located at smaller angular separations that probe 


tighter orbits (e.g. 2014 


2012 2012). In-development extreme AO sys- 
tems and upgrades to first-generation extreme AO systems 


will provide images of many exoplanets near the ice line 
Males et al.|2020). The Coronagraphic 
Instrument on NASA's Roman Space Telescope (Roman- 
CGI) could provide the first detections of mature planets 
in reflected light (Figure[2). 

Planned ground-based extremely large telescopes and 
proposed space missions promise to make the discovery and 
confirmation of a habitable, Earth-like exoplanet around a 
Sun-like star a reality within the next 25 years (e.g. 


Morales et al.|2019| (Gaudi et al.\2020}|The LUVOIR Team 
2019 2021). These facilities will endeavor to 


reveal biomarkers — e.g. water, oxygen, ozone — in individ- 
ual systems. Their surveys will provide the first assessment 
of potential habitability around stars of different masses and 
thus the true context for life on Earth. Technological in- 
novation, atmospheric characterization, and demographic 
studies of young jovian planets over the past decade pro- 
vide first, key steps towards this goal. 

In this Chapter, we provide an updated description of 
the state of our knowledge about detecting and characteriz- 
ing exoplanets by direct imaging. The last dedicated direct 
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imaging review chapter in Protostars and Planets was writ- 


ten in 2007 (Beuzit et al.|2007), before the first incontro- 
vertible exoplanet imaging detections: superjovian planets 
around HR 8799 and f Pic 
[Lagrange et al.|2010). Previous reviews from 
and [Bowler] (2016) bracket the start 


and end of the era of direct imaging surveys with facil- 
ity, conventional AO systems. Now, >5 years later, the 
first generation of extreme AO surveys have ended, bring- 
ing new discoveries, fundamentally new insights into atmo- 
spheric properties of individual systems and imaged planet 
demographics, motivating substantially more powerful in- 
strumentation and tangible plans for directly detecting true 
solar system analogues, including Earths. 

This Chapter is pedagogical in nature and is organized 
as follows. In Section 2, we outline key direct imaging 
instrumentation and methods, describing the challenges in 
achieving planet detections by direct imaging, and detailing 
critical, novel hardware and software needed to image plan- 
ets. Section 3 summarizes our current inventory of directly 
imaged exoplanets, discusses challenges with interpreting 
direct imaging detections, and outlines synergies with other 
detection techniques. Section |4| focuses on atmospheric 
characterization of exoplanets by direct imaging and spec- 
troscopy, using other substellar objects (i.e. brown dwarfs) 
as anchors to provide empirical constraints on exoplanet at- 
mospheres and employing atmospheric models to infer in- 
trinsic properties, including clouds, chemistry, and gravity. 
Next we overview the architecture of directly imaged plane- 
tary systems gleaned from astrometric monitoring, dynam- 
ical models, and planet-disk interactions (Section 5). Sec- 
tion 6 summarizes recent direct imaging surveys, the oc- 
currence rates and demographics derived from them, and 
what these results mean for models of jovian planet for- 
mation. Finally, in Section 7 we forecast the future of di- 
rect imaging, describing how technological innovations and 
new, vastly more powerful facilities may provide humanity 
with a glimpse of a habitable world for the first time. 


2. Direct Imaging Instrumentation and Methods 


Direct imaging detections require separating the halo of 
bright, highly structured scattered starlight from faint ex- 
oplanet light. The key metric used to determine the de- 
tectability of extrasolar planets is the planet-to-star contrast 
ratio at a given off-axis angular separation: the contrast 
ratio needed for a detection varies with planet properties 
(Sect 2.1). Critical hardware like AO — or (more generally) 
wavefront control systems consisting of sensors and de- 
formable mirrors — and coronagraphs sharpen and then sup- 
press scattered starlight (Sect 2.2) to achieve deep raw con- 
trasts. Novel observing techniques allow post-processing 
algorithms to further remove residual starlight, increasing 
achievable contrast ratios and thus improving planet detec- 
tion capabilities (2.3). 
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2.1. Detectability of Planets by Direct Imaging 


Even in the absence of an atmosphere to blur starlight, 
an astronomical image of a point source (e.g. a distant star) 
will not be a single point: its intensity distribution follows 
from the Fourier transform due to the telescope pupil func- 
tion. For a simple unobscured circular aperture, the image 
intensity follows an Airy pattern whose values compared to 
the peak intensity (Jo) at angular separation p, effective tele- 
scope aperture diameter D, and wavenumber k = 27/) is 
I= [goes ga , Where J; is the Bessel func- 
tion of the first kind of order one. The first null of the Airy 
function occurs at ~ 1.22 A/D,a. The full-width at half- 
maximum (FWHM) of this point source on the sky, defining 
the telescope diffraction limit, is given by 


ba) 


Even current telescopes in principal have the angular reso- 
lution needed to image planets on solar system-like scales 
orbiting nearby stars. An Earth twin orbiting at a pro- 
jected separation of 1 au from a Sun-like star at 10 pc sub- 
tends an angle of 01. Diffraction limited imaging with the 
Hubble Space Telescope at optical wavelengths or with 8- 
10m class telescopes like VLT, Subaru, and Keck at near- 
infrared (IR) wavelengths is sufficient to resolve objects at 
Earth-Sun projected physical separations out to a distance 
of 10-20 pc. However, the star's halo light must be reduced 
through hardware and software to a level sufficient to make 
the planet's light detectable. The starlight suppression level 
required to detect a planet by direct imaging depends on 
whether the detection is in reflected light or thermal emis- 
sion and varies with the planet angular separation, age, size, 
temperature, and other properties (see Figurep). 
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Fig. 2.— Contrast vs. angular separation for current and future 
ground and space-based high-contrast imaging platforms vs. mea- 
sured and predicted contrasts for selected directly imaged planets, 
radial-velocity (RV) detected planets, and Earth and Jupiter ana- 
logues. Generated using tools developed by Vanessa Bailey, which 
are available at: https://github.com/nasavbailey/DI-flux-ratio-plot. 
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2.1.1. 


Planets reflect the light of the stars they orbit. At optical 
wavelengths, the contrast ratio for a planet in reflected light 
can be approximated as 


Reflected Light 


2 
Coptical,d dd Ag (A) $(A, a) (2) ; (2) 
p 
where ry is the planet radius, ap is the planet-to-star physi- 
cal separation, A, is the visible geometric albedo spectrum, 
$ is the phase function as a function of a, the phase angle, 
which is the angle between the star, planet, and observer 
(Greco and Burrows|2015} [Traub and Oppenheimer[2010). 
A, and ¢ depend on the planet's atmospheric properties. 
For a Lambertian phase function, valid for high albedo at- 
mospheres, the phase function follows a simple relation of 
ola) = [sin(a)+(7-a)cos(@)|/7, where ¢ = 1/7 at maximum 
elongation (o = 7/2). Other potential scattering phase func- 
tions include isotropic scattering and Rayleigh scattering. 
The value of ¢ is considerably lower for Rayleigh scatter- 
ing appropriate for jovian atmospheres at angles of 60—90? 
(Madhusudhan and Burrows|2012). 
For an exo-Jupiter and exo-Earth emitting as Lamber- 
tian spheres with measured geometric albedos of 0.52 and 
0.367, these contrasts at maximum elongation reduce to: 


2 2 
5.2 
Coptical,J ~ 1.4 x 107° (2) ( =) O) 
Tj Ap 


? /1au\? 
Coptical, am 2.1 x 10-19 (2 3 (m , (4) 
Tg ap 


where ry and rg are the radii of Jupiter and Earth re- 
spectively. At V band, a Jupiter and Earth at 10 pc will 
then have apparent magnitudes of slightly greater than 27 
and 29, respectively, and require starlight removal at angu- 
lar separations of 0/5 and 0/1 better than one part in 107° 
and 10-9, respectively. These required contrasts are well 
beyond the capabilities of current ground and space-based 
high-contrast imaging instruments. Thus far, exoplanet di- 
rect imaging has therefore focused on the detection of ther- 
mal emission from self-luminous jovian planets. 


and 


2.1.2. Thermal Emission 


Like brown dwarfs, jovian planets cool and contract with 
time, releasing gravitational potential energy as thermal 
emission (e.g. Burrows et al.[2001). At ages of 1-10 Myr, 
models for the luminosity evolution of 1-10 M; exoplan- 
ets predict temperatures of ~ 500-3000 K and radii up of 
Ty ~ 2.5-3 Ry. By 1 Gyr, these models predict that these 
planets cool to temperatures T < 500 K and contract to 
Jupiter-like radii (rp ~1-1.2 Rj)(e.g. Fig. 5 in [Spiegel] 
land Burrows|2012). Thus, thermal emission from young 
jovian exoplanets peaks at near-to-mid IR wavelengths at 
the youngest ages, moving to well into the mid-IR at older 
ages (Figure Bp. The contrast ratio for planets in thermal 
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emission depends on the radius and effective temperature 
of the planet and star and other properties: 


E (5) 


where £,, (T, X) is the thermal flux from the planet as a 
function of wavelength, P .(T) is the thermal flux from 
the star as a function of wavelength, and X depends on the 
planet's atmospheric characteristics, such as clouds, chem- 
istry, and gravity. 
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Fig. 3.— Contrast ratio vs. selected optical, near IR, and mid-IR 
passbands for models for known directly imaged planets (8 Pic b 
and HR 8799 cde), models for older cooler exoplanets, and Jupiter 
assuming a G2V primary in all cases. The actual planet-to-star 
contrast ratios for 8 Pic b and HR 8799 cde are 5-10 times higher 
as they orbit intrinsically brighter, hotter stars. Plot adapted and 


modified from (2014a). Models drawn from|Currie| 
jet aL. (2011); [Madhusudhan et al.\(2011);|Burrows et al.|(2003). 


Fig. |3| shows how the contrast ratios vary for various 
types of exoplanets orbiting a G2V star. In the near-IR, 
young superjovian (~10 Mj) exoplanets similar to 8 Pic 
b and HR 8799 cde are roughly a factor of 107? to 1075 
fainter than a Sun-like star. In the red optical, predicted 
contrasts from these planets are 100 times larger; however, 
mid-IR contrasts are a factor of 10 smaller. Jovian exoplan- 
ets of similar masses around much older stars are far cooler 
— another factor of 100—1000 fainter in the red optical and 
near-IR — but remain bright in the mid IR. A true Jupiter 
analogue emits negligible thermal emission shortward of A 
74—5 um. Rocky terrestrial planets also emit thermally. 
An exo-Earth analogue with Tg œ~ 260 K has blackbody 
emission peaking at 10 jum: its contrast ratio is ~1077, a 
factor of 221000 shallower than its reflected light contrast. 


2.2. Wavefront Control and Coronagraphy 
2.2.1. Atmospheric Wavefront Control for Direct Imaging 


Wavefront errors, whether induced by atmospheric tur- 
bulence on the ground or intrinsic to an optical imaging 
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system on the ground or in space, substantially limit an im- 
age system's achievable contrast. Left unmitigated, these 
errors preclude the direct detection of exoplanets. In the 
past decade, the direct imaging field has made key strides 
in wavefront control hardware and software to drastically 
reduce wavefront errors and has developed sophisticated 
coronagraph designs to further suppress diffracted starlight. 

For ground-based imaging systems, turbulence arising 
from many atmospheric layers along the path from the 
star to the telescope induces changes in the optical path 
length difference of starlight (Guyon||2005), blurring im- 
ages. These aberrations must be corrected by an AO sys- 
tem consisting of wavefront sensors and deformable mirrors 
(DMs). An AO system splits incoming light from a guide 
star between the science detector and a wavefront sensor 
to measure and then correct — with a DM(s) — atmospheric 
turbulence distorting the incoming wavefront, sharpening 
starlight at the image plane. 

Key terms driving the wavefront error budget for an AO 
system include 1) measurement error (om) which depends 
on the noise properties of a wavefront sensor and the guide 
star brightness, 2) temporal bandwidth error (o+) which de- 
pends on the AO system time lag 7 compared to the at- 
mospheric coherence time To, and 3) the fitting error (or) 
which depends on the coherence length ro compared to the 
DM actuator density. The quadrature-added sum of these 
terms sets the Strehl ratio — a measure of the optical quality 
of the image — and the raw contrast vs. angular separa- 
tion (e.g. [Zyson]1998). At the angular separations relevant 
for most current planet searches (0 ~ 02-10), wavefront 
measurement error and temporal bandwidth error set the 
contrast floor. Wavefront chromaticity and non-common 
path aberrations due to instrument optics can also limit con- 
trast Currie et al.|2020a). 

Facility, conventional AO systems were used for the first 
direct imaging searches. They typically sample and cor- 
rect the incoming wavefront at frequencies of 0.2-1 kHz 
and use DMs with a few hundred actuators: aberrations are 
measured by standard Shack-Hartmann wavefront sensors 
(e.g. see[Tyson|1998). These systems yielded partial cor- 
rections, achieving modest Strehl ratios (SR ~ 0.1-0.4) at 
near-IR wavelengths focused on by most exoplanet imaging 
searches (e.g.|Lafreniére et al.|20077b). 

The past 5-10 years have seen the deployment of nu- 
merous extreme AO systems on 5-10m telescopes, which 
have wavefront control loop speeds of 21 kHz and DMs 
with 21000 actuators. Examples of extreme AO systems 
include the Large Binocular Telescope Adaptive Optics sys- 
tem (LBTAO), PALM-3000 on the Hale telescope at Palo- 
mar Observatory, the Gemini Planet Imager on Gemini- 
South, and the Spectro-Polarimetric High-contrast Exo- 
planet REsearch instrument (SPHERE) at the VLT, the 
Subaru Coronagraphic Extreme Adaptive Optics project 
(SCExAO) on the Subaru Telescope, and MagAO-X on the 
Clay telescope at the Magellan Observatory (Esposito et al. 


Macintosh et al.|2014 
Jovanovic et al. 015]|Males et al. 2020). These 


eral. 2019 
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systems obtain higher-quality AO corrections — Strehl ratios 
of 0.7—0.95 at 1.6 um — and yield a factor of 10-100 deeper 
contrast at sub-arcsecond separations than conventional AO 
systems (e.g Macintosh et al.|2014). 

Extreme AO systems have relied on hardware develop- 
ments in three key areas to operate high-performance AO 
correction: detectors, deformable mirrors, and computing 
hardware. Most extreme AO platforms utilize fast ultra-low 
noise detectors like Electron-Multiplying Charge-Coupled 
Devices (EMCCDs) to record stellar photons used for wave- 
front sensing with reduced measurement errors (e.g. 
Males er al 030) 
High-actuator count DMs are now available. For exam- 
ple, the microelectromechnical-type DMs (MEMS) offer 
unprecedented actuator density (7:10 act per square mm), 
allowing for fairly compact extreme AO instruments (e.g. 
[2020). Other systems employ adaptive secondary mirrors 
(usually voice-coil DMs) whose low emissivity makes them 
especially well suited for exoplanet imaging in thermal IR 
(Esposito et al.\2011). High performance computing hard- 
ware with low latency can handle the demanding computa- 
tion requirements of fast, high order extreme AO systems. 

Extreme AO systems also differ from conventional AO 
systems in their optical design and system architecture. 
Wavefront sensors (WFS) optimized for speed, sensitivity, 
accuracy and precision such as the Pyramid wavefront sen- 
sor have been adopted by many leading extreme AO sys- 
eme op panel of igre peste ero noe 
Males et al.|2020), while the more established 
Shack-Hartmann WFS approach has been upgraded with 
spatial filtering to improve performance. Wavefront correc- 
tion is often performed in two steps: a conventional coarse 
"woofer" correction is followed by a faster and more ac- 
curate “tweeter” correction (e.g. [Macintosh et al.|2014). 
Extreme AO systems also include dedicated sensing and 
control of pointing for precise co-alignment of the star and 


coronagraph mask (Singh et al.|2015||Shi et al.|2017 
[et al. 2017). which is essential to maintain high contrast, 


especially at the smallest angular separations. 

Advances in wavefront control algorithms are also en- 
abling significant improvements in high-contrast imaging 
performance. Recent advances include predictive control 
and sensor fusion. These and other new software/hardware 
advances are described in Section 7.1. 


2.2.2. Optical Starlight Suppression 


Starlight suppression is primarily performed by dedi- 
cated coronagraph optics: phase and amplitude masks de- 
ployed in the beam to remove starlight while preserving 
planet light. These optical elements include occulters in 
the image plane to directly block out on-axis starlight, as 
well as optical elements in the pupil plane to manage the 
telescope diffraction pattern. Many coronagraph designs 


and approaches are available (e.g.|Kuchner and Traub|2002 
2003| |Sivaramakrishnan et al.|2005 
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Fig. 4.— Some key high-contrast imaging hardware compo- 
nents. (top panel) Wavefront control with a Pyramid WFS and 
high-performance camera: images of the 4 pupils created by a 
pyramid optic and acquired by a photon counting, high speed, 


low-latency EMCCD 2015). (bottom panel) 


Advanced starlight suppression masks: electron scanning micro- 
scopic image of an annular groove phase mask (“vector vortex” 


coronagraph) (Mawer er al. 012] [Delacroix er al 2013). 


2006| |Kenworthy et al.12007 2010 
2017), providing a wide range of performance 


characteristics. For example, the well-established conven- 
tional Lyot coronagraph combines an occulter in the image 
plane with a Lyot stop in the pupil plane to deliver robust 
performance at moderate contrast and large angular separa- 
tions. Phase-mask coronagraphs or interferometric designs 
(e.g. the vector vortex coronagraph) can be highly opti- 
mized for deep contrast at small inner working angles (Fig- 
ure[4] bottom panel). The highest-performing coronagraphs 
are also the most demanding in terms of wavefront qual- 
ity and stability. Consequently, extreme AO systems first 
adopted Lyot coronagraphs, while more recent systems or 
upgrades deploy higher performance solutions as the corre- 
sponding wavefront quality improves. 

Wavefront control is also becoming an important com- 
ponent of starlight suppression in the form of speckle con- 
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trol. In speckle control, a feedback loop from the post- 
coronagraph image to an upstream deformable mirror al- 
lows for residual speckles to be measured and canceled. 
This approach has been used in laboratory testbed to reach 
deep contrast levels (1077—107!9) in support of future 
space-based missions. Deployment on ground-based sys- 
tems remains challenging due to the dynamic nature of at- 
mospheric turbulence, but has been successful in removing 
a fraction of static and slow speckles (Codona and Kenwor-! 


let al.||2021). Sensing and compensation of non-common 
path aberrations can also be performed using a dedicated 
sensor, e.g. a Zernike phase-mask sensor 
2019}. 


2.3. Observing, Post-Processing, and Spectral Extrac- 
tion Methods 


At sub-arcsecond separations, the raw (instrument- 
delivered) planet-to-star contrasts obtained by the first gen- 
eration of facility and now extreme AO systems — 107 3— 
1074-5 — are still too shallow to yield decisive detections 
of many young jovian exoplanets. Residual quasi-static 
speckle noise, due to nanometer-scale imperfections in 
telescope optics, thermal flexure, etc., limits instrumen- 
tal contrasts and fails to simply average out over time like 
photon (“white”) noise (Hinkley et al.[2007). Speckle noise 
statistics follow a modified Rician distribution whose long 
positive tail can lead to many false positives 
[Soummer et al.|2007). Therefore, novel observing 


techniques coupled with advanced post-processing methods 
are critical to optimizing our ability to image exoplanets by 
further suppressing stellar halo light by orders of magnitude 
and making residual noise more Gaussian (“whitening”). 

Various forms of differential imaging techniques are uti- 
lized to remove speckle noise for direct imaging. Com- 
monly used techniques include: 


* Angular Differential Imaging (ADI) 
2006) — ADI exploits the quasi-static nature of resid- 


ual speckle noise in AO corrected high-contrast im- 
ages from the ground and diffraction-limited imaging 
from space. For ground-based AO imaging, quasi- 
static speckles evolve on a characteristic timescale 
of ~ 10s of minutes to an hour or more 
et al.\2007). For space-based imaging, PSF breathing 
due to thermal variations and other small mechani- 
cal changes cause the quasi-static halo to evolve, but 
the halo for a given target can remain well corre- 
lated for months or years. By turning off the instru- 
ment rotator that keeps north fixed on the detector (on 
the ground) or obtaining observations at different roll 
angles (space), astrophysical objects off-axis from a 
star change position angle on the detector while the 
speckle halo remains fixed. Subtracting images ob- 
tained in an ADI sequence can then suppress speckle 
noise without fully suppressing planet signals. 


ADI efficiently suppresses speckles on the ground 
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and in space (Marois et al.|2006). It is the most 


widely used observing method in direct imaging. 
ADI’s advantages are most limited at small angular 
separations where the displacement of the planet's 
PSF is smallest for a given parallactic angle change or 
roll angle change. Typically, ADI induces signal loss, 
which must be corrected to achieve absolute spec- 
trophotometric calibration for planets (see below). 


Spectral Differential Imaging (SDI) (Marois et al.| 
Biller eral 204) — At 
raw contrasts shallower than the 107‘ levels needed 
to image self-luminous jovian exoplanets, achromatic 
phase errors dominate the wavefront errors respon- 


sible for the speckle halo (e.g. |Shaklan and Green 
2006 2018). For data obtained simulta- 


neously at multiple wavelengths — e.g. dual-channel 
imaging or an integral field spectrograph (IFS) - the 
speckle halo as a function of wavelength, scaled in ra- 
dius by the central passband wavelength, is extremely 
well correlated. Since the images at different wave- 
lengths are obtained simultaneously with the target, 
they do not suffer from are temporal decorrelation as 
with ADI (e.g. [Marois et al.|2014). SDI then sup- 
presses the speckle halo in a given passband by con- 
structing a reference PSF drawn from (rescaled) im- 
ages at other wavelengths. 


Like ADI, SDI induces signal loss and is least ef- 
fective at small angular separations, where band- 
pass rescaling yields small displacements of a (de- 
)magnified planet PSF. SDI's efficacy is also limited 
by non-common path errors, differential sampling of 
the point-source PSFs across wavelength channels, 
and the spectrum of the planet to be detected (e.g. 
ing throughput-corrected spectra from data processed 
with SDIis also more challenging than with other dif- 


ferential imaging techniques (Pueyo|2016). 


Reference Star Differential Imaging (RDI) — RDI 
subtracts the PSF of a companion-less reference star 
from the target. It is widely used in the ultra-stable 
environment of space. While ADI and SDI may yield 
superior speckle suppression at moderate to wide sep- 
arations, RDI may be advantageous at very small an- 
gular separations where other techniques suffer from 
self subtraction (Xuan et al.|[2018). RDI requires a 
very stable PSF to efficiently operate. Using it on 
ground-based telescopes may require fast switching 
between the target and reference star(s) for maximum 
contrast gain (e. g.[Wahhaj et al.[2021); in broad band- 
passes, it also requires a reference star that is ex- 
tremely well color matched (e.g. |Krist et al.|1998). 


These differential imaging techniques are used in combi- 
nation with PSF subtraction algorithms. Many widely-used, 
advanced PSF subtraction methods are different forms of 
least-squares algorithms, which construct a reference PSF 
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No AO correction Extreme AO correction 


that minimizes the variance when subtracted from a given 
target image. The locally-optimized combination of images 
(LOCD algorithm constructs a reference PSF from a linear 
combination of reference images weighted by coefficients 
c determined from the solution to the matrix equation c = 
Atb (Lafreniére et al.|/2007a), where A is the (square) 
covariance matrix for the reference library and b is the col- 
umn matrix populated by elements multiplying each ref- 
erence image by the target image. Another approach ex- 
ploits principal component analysis (PCA), computing the 
Karhunen-Loéve transform of the reference image set and 
projects this set onto the target image to construct a combi- 
nation of weighted eigenimages to subtract from the target 
(Karhunen-Loéve Image Plane algorithm, henceforth KLIP, 
Amara and Quanz[2012). Multi- 
ple successor algorithms such as TLOCI and A-LOCI draw 
from the lineage of LOCI and/or KLIP, employing advances 
such as correlation-based frame selection, pixel masking, 
various rank-truncations of the covariance matrix, free pa- 
rameter optimization or exploiting high-performance com- 
puting to solve for the variance-minimizing coefficients di- 


rectly (eg. 2014} Currie ei al 0128] 
Thompson and Marois|2021). Separate approaches 
involve maximum-likelihood methods to model and remove 
the stellar PSF (e.g. Cantalloube et al.|2015) or statistically 
modeling non-stationary covariances in small regions of im- 
ages to improve PSF subtraction (e.g. the PACO algorithm 
Flasseur et al.|2018). 

The relative performances of these algorithms vary in the 
literature and may depend on their exact implementation 
and suitability for a particular data set, although some work 
suggests that the newer algorithms descended from LOCI 
and KLIP can offer significant improvements (e.g. Rameau| 
[Thompson and Marois 021). In general, all 
of these algorithms enable significant contrast gains over 
simple, classical methods, especially at small angular sepa- 

ADI, SDI, and RDI combined with PSF subtraction 
methods can attenuate and distort planet signals. The first 
methods to measure these biases and recover true planet 
flux measurements and astrometry injected synthetic plan- 
ets into data at other locations or iteratively subtracted 
negative copies of planets at their apparent positions (e.g. 


2007a| Lagrange et al.|2010). Recently, 


several authors have developed efficient forward-modeling 
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* PSF subtraction 


* coronagraph 
Fig. 5.— A demonstration of high-contrast imaging hardware and software with SPHERE. Seeing limited point-spread function (PSF) 
(left) sharpened by an extreme AO system (middle left) with scattered starlight further suppressed by a coronagraph (middle right) and 
advanced PSF subtraction methods (right). 


methods to estimate photometric and astrometric biases at 


the planet's location (e.g. 2016). Planet forward- 


models in turn can be used as matched filters to improve 
planet detection capabilities themselves (Rufio et al.|2017). 

Analyses of images whose speckle noise is partially sup- 
pressed by advanced PSF subtraction techniques have re- 
vised notions of how to quantify detection significances 
and spectroscopic uncertainties. For instance, finite sample 
sizes impact our definitions for contrast limits drawn from 
images with noise whitened by KLIP, LOCI, and other al- 
gorithms, especially at small (1-3 A/D) separations 
[et al.|2014} |Pairet et al.[2019). In IFS data, residual speckle 


noise may be spatially and spectrally correlated. Consider- 
ing the full spectral covariance has a substantial impact on 
deriving planet atmospheric parameters from model com- 


parisons (Greco and Brandt|2016). 
Figure |5| shows the combined effect of high-contrast 


imaging hardware and software on imaging exoplanets, 
stepping through the successive improvements found from 
using extreme AO, coronagraphy, and PSF subtraction. 


3. Direct Imaging Detections 


3.1. 


Figure |1| shows the demographics of planetary mass 
companions to stars detected via various methods. Commu- 
nity consensus on the planethood of many imaged objects 
(e.g. HR 8799 bcde, 8 Pic bc, 51 Eri b, etc.) is clear. How- 
ever, identifying the exact criteria needed to distinguish be- 
tween planets and brown dwarfs is challenging. 

Planets have often been identified as objects with masses 
below the deuterium-burning limit, nominally 13 Mj (e.g. 
the IAU Working Group Definition). However, this sim- 
ple criterion is poorly motivated. In addition to being time 


and metallicity dependent (Spiegel et al.|2011), deuterium 


burning arguably does not identify a meaningful bound- 
ary for the evolution of low-mass objects 
Luhman|2008), as some objects below the deuterium 
burning limit have been found in configurations which im- 
ply formation by cloud formation and other objects above 
the deuterium burning limit have been found in configura- 
tions which imply formation like a planet in a disk. For in- 
stance, recent imaging surveys have found objects that are 
members of quadruple systems, clearly formed by molec- 
ular cloud fragmentation, with inferred masses down to 5 
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M; and free-floating objects with sub-deuterium burning 
masses as well [Liu et a.[2013). In 
contrast, RV surveys have identified some systems — e.g. 
the 2.7 Mo star v Oph — with companions at ~ 1 au with 
masses of 22 and 24 Mj that are nevertheless locked in a 
mean-motion resonance indicating formation in a disk (i.e. 
like a planet) (Quirrenbach et al.|2019). The physics of 
planet formation does not require that gas accretion shuts 
off once 13 M; of material is accreted. 

An alternate definition leverages formation processes: a 
planet is an object formed in a circumstellar disk around a 
young star. Demographic analyses of substellar objects can 
provide empirically-motivated criteria for separating plan- 
ets from brown dwarfs. Studies of the substellar mass func- 
tion from previous RV surveys and recent ones (i.e. the Cal- 
ifornia Legacy Survey) show a local minimum at msin(i) ~ 
16-30 My (e. 
T. Currie et al. submitted). The minimum in the com- 
panion mass function may be proportional to the primary 
mass, indicating that companion mass ratio (q) could be a 
key discriminator (Grether and Lineweaver|2006). Theory 
also suggests that the semimajor axes (ay) and mass ratios 
(q) of companions also help distinguish between bona fide 
planets and brown dwarfs 
2011). Binary companions to more massive stars with q = 
0.025 are exceptionally rare 
[et al.|2016). Distributions of protoplanetary disk radii peak 
at ~200 au and fall to low frequencies by ~300 au (e.g. 
[Andrews and Williams|2007): companions at wider separa- 
tions are far less likely to have formed from a disk unless 
scattered to their current locations by unseen companions. 
Thus, we set the following limits for a planet vs. a brown 
dwarf: mass < 25 Mj, q < 0.025, and ap 5 300 ad 

Our linked spreadsheet lists the current inventory of 
directly imaged exoplanets, protoplanets, higher mass- 
ratio/wider separation planet-mass companions, and con- 
troversial caseg] As of 1 May 2022, we identify 20 directly 
imaged exoplanets that fit our criteria, three of which are 
imaged protoplanets. We also list another 34 dozen com- 
panions that may instead be better considered as brown 
dwarfs instead of planets (including companions orbiting 
brown dwarfs), and 6 “controversial” cases. Some of these 
individual classifications will undoubtedly change pend- 
ing new analysis as may the exact values used to separate 
planets from brown dwarfs. Objects in each category will 
certainly be added over the next few years. 

Below, we describe general properties of fully-formed 
directly imaged planets within 300 au and discuss three 
well-studied, emblematic cases (HR 8799, 8 Pic, and 51 
Eri). We also summarize our current knowledge of proto- 
planets and challenges with identifying imaged exoplanets. 


1Our mass limit is similar to that adopted elsewhere: e.g. the NASA Exo- 
planet Archive. The mass ratio and separation criteria further restrict our 


sample. 
^https://tinyurl.com/srb33b 
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32. Fully-formed Exoplanets 


3.2.1. General Properties 


Most directly imaged planets known thus far have near- 
IR contrasts ranging between 1074 (e.g. xk And b) and 107° 
(51 Eri b). All of them are self-luminous, made visible by 
thermal radiation released as the planets cool and contract. 
On the sky, the planets generally lie at angular separations 
of 072 to 2/0 (~ 5-40 A/D for an 8m telescope). 

Stars hosting imaged exoplanets typically have ages of 
10-100 Myr; many are found in young nearby stellar asso- 
ciations (often «100 pc, with ages «100 Myr) that share 
similar kinematics (proper motions, space motions) and 
formed in the same star-forming region, such as Sco-Cen 
or the 8 Pic Moving group (e.g. 
Song et al [2003] [Schlieder er al. 
! 2012). (see 


Section [6.1] for more details on such associations as sites 
for exoplanet searches). Most stars with imaged planets to 
date are B, A, or F stars, at least 5096 more massive than 
the Sun. Most systems with imaged planets — including sys- 
tems with the first detections (HR 8799, 8 Pic) — also have 
Kuiper belt-like debris disks. 

Planets detected to date share similar spectral types and 
temperatures with more massive field brown dwarfs, but 
often show evidence of low surface gravity in their spec- 
tra (Barman et al.[2011b[a). Brown dwarfs and exoplanets 
alike cool monotonically with age, beginning life as hot M- 
type objects, cooling to the L spectral type (with very red 
near-IR colors and silicate condensate clouds 
2005), the cooler T spectral type (with blue near-IR col- 
ors and strong methane absorption at 1.6 wm and 2.2 um, 
Kirkpatrick|2005) and eventually to the very cool Y spectral 
type (Cushing et al.[2011} [Kirkpatrick et aL[2012). Thus, 
there is an age / mass / temperature degeneracy for these ob- 
jects, rendering mass estimates based on models very sen- 
sitive to the age assumed for the system. 

Most mass estimates for imaged planets depend on lumi- 
nosity evolution models, challenges for which are described 
in Section 3.4] Some planets — e.g. HR 8799 bcde, 6 Pic b, 
51 Eri b, and HD 206893 b — have dynamical mass measure- 
ments or limits. Masses inferred from luminosity evolution 
or derived from dynamics are typically ~5-10 My. Some 
planets have masses near the deuterium-burning limit (& 
And b, HD 206893 b); 51 Eri b may have a far lower mass 
(as little as ~2 Mj). Mass ratios for most imaged planets 
are q ~ 0.005-0.01: 51 Eri b may have the lowest mass ra- 
tio (q ~ 0.001). Even leading extreme AO systems are typi- 
cally not sufficiently sensitive to detect young Jupiter-mass 
planets, let alone Saturn-mass planets. However, future ca- 
pabilities will close these gaps (See Section[7). 


3.2.2. Emblematic Systems 


While planetary mass companions to brown dwarfs 
had previously been detected as early as 2004 (most no- 


tably 2M1207b, 2004), the near simultane- 


ous announcements of HR 8799bcd and £ Pic b are widely 
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Fig. 6.— Well-studied directly-imaged exoplanetary systems: (left) HR 8799, (middle) 8 Pic, and (right) 51 Eri. Images draw from 


(2010b), |Lagrange et al.| (2019), and (2015), respectively. Note that 8 Pic c is not visible in the middle 


panel: the HR 8799 and 51 Eri images show the full complement of known planets in each system. 


regarded as the first bonafide directly-imaged exoplanet 
detections [Lagrange et aL 2009) 
Planets in both systems were detected from ground-based 
facility AO systems in the near-to-mid IR in thermal emis- 


sion, not scattered light?| Just over 7 years later, Macintosh 
g y 


(2015) announced the first exoplanet discovered with 
extreme AO: 51 Eri b. 


HR 8799 — HR 8799 is a nearby (~40 pc) mid-A field 
star with an estimated age of ~40 Myr (Baines et al.[2012). 
The star hosts a massive, resolved Kuiper belt-like debris 
disk and a warm debris population interior to 10 au con- 


sistent with an asteroid belt (Sadakane and Nishida|1986 
[Su et al.[2009 2014] 2016): 


possible signposts of massive, perturbing planets. 


In 2008, |Marois et al.| announced the direct imaging dis- 
covery of HR 8799 bcd followed by a fourth planet dis- 


covery in 2010 (HR 8799 e; [Marois et al.|2010b), located 
at projected separations of 15-70 au (Figure 6| left panel). 
Soon after HR 8799 bcde's announcements, other studies 
identified one or more planets in archival or separately- 
obtained data 
[2012a). A decade of Keck Observatory monitoring 
showed that HR 8799 bcde orbit close to the 1:2:4:8 res- 


onance (Konopacky et al.12016). Current analyses suggest 


that the planet orbits are nearly coplanar with the disk, with 


a small inclination of ~27° (Konopacky et al.12016). How- 


ever, HR 8799 e may not orbit on a plane strictly coplanar 
with HR 8799 bcd (Gravity Collaboration et al.|2019). 
Originally, masses inferred from the planets' luminosi- 
ties spanned a wide range of values (5-13 My) 
[2010a). However, dynamical stability mod- 
eling strongly favors masses below 10 M; (7 My) for HR 
8799 cde (HR 8799 b) 


2011} |Sudol and Haghighipour|2012). Using Hipparcos 
and Gaia, (2021b) measured a dynamical mass 


for HR 8799 e of 9.6* '$ Mj, consistent with these lim- 
its (see Section B.5). Additional planets may explain the 


2 Another object — Fomalhaut b 2008) — was announced on 


the same day as the HR 8799 planets but appears to be made visible by 
scattered starlight (Currie et al.|2012b||Galicher et al.|2013). Its nature is 
further discussed in Section 


shape of the inner edge of HR 8799's cold belt (Booth et al.| 


2016), especially if they are near or below Jupiter's mass 
TORIS LE MIB ARR Any planets interior 
to HR 8799 e must be below 3-4 M; at 7-10 au and 5-6 
M; at 4-7 au Brandt et al.|2021b). The 


HR 8799 planetary system resembles a scaled-up version of 
our own outer solar system (Marois et al.[2010a). 

HR 8799 bcde have been benchmark objects for under- 
standing the atmospheres of young jovian planets. Their 
photometry and low-resolution spectra differ from those of 
older, field substellar objects thought to have similar tem- 
peratures, identifying features diagnostic of clouds, chem- 
istry, and gravity (Currie er aL [2011| 
2014b). Higher resolution spectra probe molecular abun- 
dances connected to formation mechanisms 
[Barman et al.[2015). Section 4] discusses HR 


8799 bcde's spectra and atmospheres in more detail. 

B Pic — The A6V star 8 Pictoris is the eponymous mem- 
ber of a collection of kinematically-associated stars known 
as the 8 Pic Moving Group (age~20 Myr). The star hosts 
an edge-on debris disk, imaged since the mid-eighties from 
the ground and space (e.g. [Smith and Terrile[1984). The 
disk consists of a planetesimal ring at 70—100 au responsi- 
ble for the bulk of the dust through collisions and an outer 
extension comprised of small dust grains blown out by ra- 
diation pressure (Lagrange et al.|2000). 

enr TODOS Rcs knife B Pic b from data 
obtained in 2003 with VLT/NaCo; follow-up data obtained 
one year later confirmed the planet (Lagrange et al.|2010). 
Thanks to 8 Pic b's small semi-major axis (~ 9 au), over 
7596 of its orbit has been monitored 


Nowak et al.|2020). More recently, a second planet (8 Pic c) 
was identified from RV data and then recovered with high- 


contrast inteferometric imaging with VLTI/GRAVITY 
(Gravity Collaboration ei al. P020) 

DB Pic b may partially explain the observed warp in the 
edge-on debris disk and the many evaporating exocomets 
identified over the past 30 years 
Kiefer et al.|2014). Planets besides 9 Pic bc could explain 
why the inner 70 au is relatively devoid of debris dust; lo- 
calized features in the disk also suggest additional planets 
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0.25"/28 au 


ALMA 
(submm) 


could be present. From RV and astrometric data, 
ler al. (20214) estimate a mass of 9.3 pos My and 8.3 + 1 
M; for B Pic b and c, respectively (see also T, 
Section [3.5), although 
these values heavily depend on the assumed uncertainties 
on the RV data. RV and direct imaging data combined ex- 
clude additional planets more massive than 2.5 M jup from 
0.1 au to hundreds of au (Lagrange et al.|2020). 8 Pic b is 
more luminous and hotter than the HR 8799 planets (Teg ~ 
1700-1800 K), just slightly redder than field early L-type 
dwarfs, and is likely cloudy/dusty with a low gravity (Bon- 


nefoy et al.|2013| |Currie et al.|2013; |Chilcote et al.|2017); 


B Pic c probably has a temperature intermediate between its 
sibling and the HR 8799 planets (Nowak et al.|2020). 

51 Eri — 51 Eridani is a 29 pc-distant early F star, a mem- 
ber of the ~20 Myr-old 8 Pic Moving Group, and a mem- 
ber of a wide hierarchical triple system that includes an M- 


dwarf binary GJ 3305 (Feigelson et al.12006). 51 Eri has a 


detected infrared excess, modeled with a cold dust belt lo- 
cated approximately between 5 and 80 au 
Riviere-Marichalar et al.)2014). 

The GPI campaign (GPIES) team discovered 51 Eri b, 
a faint planet at rproj ~ 13 au (Figure 6, right panel; Mac-| 
[intosh et al.|2015). The planet is the first discovered using 
extreme AO and the first incontrovertible T dwarf planet, 
showing strong methane absorption in H band. The planet 
is likely more eccentric than either 8 Pic b or HR 8799 bcde 
e ~ 04-0,62). If 
confirmed, an eccentric orbit could indicate the presence of 
an additional massive body or could be due to gravitational 
perturbations from GJ 3305 AB. Assuming a hot-start lumi- 
nosity evolution, current data rule out other planets more 
massive than 4 M; beyond 5 au and more massive than 2 


My; beyond 9 au (Samland et al.\2017). 


VLT/MUSE 
(Ha, optical) 
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Fig. 7.— (left) ALMA image of PDS 70 showing a cleared protoplanetary disk cavity and a circumplanetary disk around PDS 70 c and 
H-a image showing the detection of PDS 70 b and c (Benisty et al.|2021 2019). (right) SCExAO/CHARIS image (rescaled 


by stellocentric distance) of AB Aurigae showing a bright, concentrated emission source revealing the location of a protoplanet 


jet al.|2022). 


The mass of 51 Eri b is not well constrained: values de- 
rived from comparing 51 Eri b’s luminosity and age to evo- 
lutionary models favor ~2 Mj, while atmospheric model- 
ing may favor larger values, up to ~ 9 My 
Samland et al.|[2017). Different characterization stud- 
ies also find slightly diverging atmosphere properties, il- 
lustrating the challenge associated with characterizing very 
faint exoplanets with direct imaging 
[Rajan et al.|2017). Analysis of Gaia and Hipparcos astrom- 


etry set an upper limit of 11 Mj for 51 Eri b (Dupuy et al. 
2022). The planet likely has a temperature of ~700 K and 


either lacks clouds or is only partially covered by clouds 


(Rajan et al./2017\\Samland et al.|2017). 


3.3. Protoplanets 


Direct images of planets in active assembly (protoplan- 
ets) around stars that still retain gas and dust-rich proto- 
planetary disks clarify how and where planets form. The 
large distances to the nearest star-forming regions (7150 
pc) mean that protoplanets orbiting their host stars at solar 
system scales are located at very small angular separations. 
However, protoplanets can be bright (L ~ 10 ?—10 ^? Lo), 
especially if they are surrounded by their own circumplan- 
etary disks (e.g.|Zhu[2015]. 

The 5 Myr-old 0.87 Mo star PDS 70 hosts the first in- 
controvertible detections of jovian protoplanets: PDS 70 b 
and PDS 70 c (Keppler eral. 2018 
Haffert et al.12019) (Figure|7| left). Both protoplanets are 


located within the PDS 70 disk cavity, at angular separa- 
tions of p ~ 0718 and 0724 and estimated semimajor axes 
7-20 and 34 au (Wang et al.|2021). Dynamical arguments 
strongly favor masses less than 10 M; for PDS 70 b, while 
PDS 70 c's mass is more poorly constrained. Masses in- 
ferred from SED modeling range between 1 and a few jo- 
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vian masses (Stolker et al.12020). PDS 70 b is slightly ec- 
centric (e ~ 0.17 + 0.06), while PDS 70 c's orbit is consis- 


tent with being circular. The protoplanets' IR data are best 
fit by model atmospheres with substantial dust/extinction. 
PDS 70 bc show Ha emission consistent with accretion 
at rates of 1-2x 1075 Mo yr, slightly less than the stellar 
accretion rate (Haffert et al.|2019), but lack evidence for 
Br-y n PDS 70 c shows direct 
evidence for a circumplanetary disk with an estimated mass 
of 0.007—0.03 Mẹ. Thermal IR data may suggest that PDS 


70 b is surrounded by a circumplanetary disk (Stolker et al. 


2020}|Wang et al.|2021). 
Recently, data from Subaru/SCExAO and the Hubble 


Space Telescope over 13 years reveal evidence for a wide- 
separation (~93 au) embedded protoplanet around the 1— 
3 Myr-old, 2.4 Mo star AB Aurigae (Currie et al.|2022). 
AB Aur b is consistent with a protoplanet responsible for 
the millimeter dust cavity and CO gas spirals both seen by 
ALMA (Tang et al.|[2017). It appears spatially extended, 
plausibly due to light from the central source reprocessed by 
the star's protoplanetary disk. The best-fit composite model 
explaining AB Aur b's optical to near-IR emission includes 
a 2.75 Ry, 9 M; source with a surface gravity of log(g) = 
3.5, emitting at a much hotter temperature than PDS 70 bc 
(72000—2500 K), and accreting at a rate of M ~1.1x107ê 
My yr- 1. The source is detected in Ha although it is un- 
clear whether this emission is due to accretion. Embedded 
in a massive disk with numerous spiral arms at over three 
times Neptune's disk instead of in a fully cleared cavity like 
PDS 70 bc, the properties of AB Aur b may point to a planet 
formation mechanism by disk instability (see Section 6). 

Prior to the discovery of PDS 70 bc, other studies 
claimed detections of protoplanets located within the gaps 
of or embedded in disks around young stars. The 2-solar 
mass, protoplanetary disk-hosting star HD 100546 has a 
protoplanet candidate at ~ 50 au and another at ~ 13 
au, just interior to the gap in the protoplanetary disk (HD 
100546 be: [Currie et aL]]OIS). HD 
100546 b has been detected in multiple data sets but ev- 
idence for orbital motion is not yet clear 
[Sissa et a.|2018); HD 100546 c has been imaged 
in a single data set and inferred through spectroastrometry 
but not yet imaged in subsequent data 
[Sissa er af 2018). Interpreting both 
candidates — whether a planet or disk feature — is challeng- 
ing (Curie et al [2015 
et al.\2018). Both candidates require further study and con- 
firmation. 

Two studies presented detections of protoplanets around 
the young, Sun-like star LkCa 15 through a combination of 
sparse aperture masking interferometry (SAM) and H, dif- 


ferential imaging (LkCa 15 bcd; |Kraus and Ireland |2012 
Sallum et al.|2015). However, later direct imaging obser- 


vations showed that the SAM detections correspond to disk 


features (Currie et al. 2019): LkCa 15 b technically remains 


a candidate for now due to its single epoch H,, detection 


(though see 2018). Other claimed proto- 
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planet detections have been revealed to likely be misidenti- 


fied disk signals instead of planets (e.g. 2019). 


3.4. Challenges with Interpreting Detections 


3.4.1. Confirming Companionship 


Direct imaging observations reveal many point sources 
that are unrelated background stars instead of bound com- 
panions, especially for systems in the Galactic plane (e.g. 
[Galicher et al.|016). Confirming candidate planets as 
bound — i.e. sharing common proper motion with and orbit- 
ing their stars — often requires multi-year observations, de- 
pending on the star's proper motion. In some pathological 
cases, background stars can have a non-zero proper motion 
and thus are more easily confused with bona fide planets 
(e.g. Nielsen et al.|2017). In the absence of full confirma- 
tion, an object's near-IR spectrum can provide strong evi- 
dence that it is a directly-imaged exoplanet (e.g. for 51 Eri 


b; [Macintosh et al.[2015). 


3.4.2. Estimating Accurate Masses 


Masses for imaged planets and planet candidates are typ- 
ically not directly measured but are instead estimated using 
luminosity evolution models that map between an object’s 
brightness and mass as a function of age. Accurate mass es- 
timates therefore are affected by uncertainties in these mod- 
els and require precise system ages. Hot start models for 
planet luminosity evolution assume a high initial entropy, 


resulting in bright planets for the first 1-100 Myr 
2002} [Burrows et al.|2001] 2003). Cold start mod- 


els assume a low initial entropy, resulting in planets that 
are substantially fainter for the first 100 Myr 
[2007). Originally, hot start models were used to describe 
planets formed by disk instability, while cold start models 
described planets formed by core accretion 
However, recent models show that core accretion- 
formed planets can be compatible with hot start-like lumi- 
nosity evolutions, and jovian planets may form with a range 


of initial entropies (Berardo et al.|2017| |Spiegel and Bur- 
[rows|2012). Aside from 51 Eri b, all planets imaged thus 


far are only consistent with a hot-start luminosity evolution. 

Precise ages can be exceptionally challenging to derive 
for isolated stars (e.g. Soderblom et al.12014). Members of 
young moving groups or other nearby associations can be 
age-dated using a variety of methods: e.g. fitting group 
Hertzsprung-Russell (HR) diagram positions with stellar 
evolutionary models, Lithium abundances. Thus, moving 
group or association membership can yield a star's age. 
However, even stars with a motion similar to bona fide mov- 
ing group members could instead be interlopers, especially 
if they have dissimilar space positions. For isolated Sun- 
like stars, stellar rotation and activity can give age esti- 


mates, albeit with significant scatter (Barnes|2007 
2011 2019). For early-type stars, HR di- 


4See Section|6.4]for definitions of the core accretion and disk gravitational 
instability formation pathways 
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agram positions have provided approximate ages: optical 
interferometry can now provide more precise estimates by 
resolving the stars themselves (e.g. Jones et al.|2016). 
Revisions in the stellar age often affects the interpreta- 
tion of an imaged companion. For example, based on the 
primary's claimed membership in the Columba association, 
& And b was originally thought to be 12-13 M; 


2013). In contrast, the primary's HR diagram position 
resembles an older system (220 Myr), which would imply 


the companion is ~ 50 M; (Hinkley et al.|2013). However, 
through optical MM E M showed 
that the star is likely (nearly) coeval with Columba even 
if it is not a member (—47 Myr), implying the object is a 
planet-mass companion. Later, near-IR spectra of « And 
b were found to be consistent with a planet interpretation 


(Currie et al.|2018). In contrast, GJ 504 B was originally 
announced as a 3-8.5 Mj planet with an age of 100-510 


Myr (Kuzuhara et al.|2013). However, subsequent analysis 


of the star suggested a revised age of ~2.5 Gyr and thus a 


mass into the brown dwarf regime (D 'Orazi et al.|2017). 
3.4.3. Planet or disk feature? 


Distinguishing between highly-structured disk signals 
and bona fide protoplanets is a steep, chronic challenge. 
Advanced algorithms needed to detect protoplanets can at- 
tenuate both disk and planet signals: forward-modeling is 
required to ensure that a claimed planet detection is not a 
partially subtracted piece of the disk (Currie et al.|2017). 
For systems with hot dust (T ~ 1000-2000 K) near the star 
that intercepts and reprocesses emission, the disk’s scat- 
tered light spectrum may strongly resemble spectra of bona 


fide protoplanets (Mulders et al.|2013). Orbital motion over 


multi-year timescales can better establish that a signal is 


an orbiting planet and not a static disk feature (Keppler 
et al.|2018). Ha emission can also pinpoint actively accret- 


ing protoplanets (Haffert et al.|2019); however, accretion 
onto protoplanets embedded in disks may be unidentifiable, 
since optical extinction likely renders Ha undetectable. 
Fomalhaut b may represent the first of yet another class 
of objects with a challenging interpretation 
(2008). The object was initially identified as a directly im- 
aged planet, made visible by both reflected light from a cir- 
cumplanetary disk and thermal emission, and responsible 
for sculpting the star’s Kuiper belt-like debris disk. How- 
ever, later analysis showed that Fomalhaut b’s spectrum 
is completely explained by scattered starlight and its or- 
bit likely crosses the ring (Janson et al.)2012}|Currie et al. 
ject is lower in mass and likely made visible purely by cir- 
cumplanetary dust. Recently, from analyzing archival and 
unpublished data, proposed that 
Fomalhaut b may be fading and dispersing, consistent with 
a massive planetesimal collision. Other analyses find no 
clear evidence for these two trends (Currie et al.[2012b), al- 
though they were conducted only over a subset of available 
data. Future observations of Fomalhaut b and/or reanalyses 
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of recent data may clarify its true nature. 


3.5. Synergies with Other Techniques 


The limitations of direct imaging affect our ability to in- 
terpret data for individual objects and to draw conclusions 
from large-scale population studies. For individual objects, 
mass estimates depend on both the age of the planet and the 
chosen luminosity evolutionary model (see Section[3.4). At 
a population level, the angular resolution of the telescope 
and the achievable contrast of the instrument limit the range 
of detectable planetary masses and semi-major axes with di- 
rect imaging. Combining direct imaging with other planet 
detection techniques partially mitigates these limits. 

Relative astrometric measurements of a planet around 
its host star cannot alone be used to precisely measure the 
mass of the orbiting planet; the semi-major axis of the rel- 
ative orbit encodes the system’s total mass, not the individ- 
ual components’ masses. The system’s mass ratio, and thus 
the mass of the planet, can be measured if the semi-major 
axis of the orbit of the star around the system barycenter 
is known. This orbit can be measured either using ab- 
solute astrometric measurements from catalogues such as 
Hipparcos and Gaia, or spectroscopic observations to mea- 
sure the Doppler shift of the star’s spectral lines over the 
course of the orbit. In some cases, high-resolution spec- 
troscopy can measure the Doppler shift of the planet’s spec- 
tral lines (Snellen et al.[2014] Wang et al.|2021), also yield- 
ing the mass ratio. Combining the star’s orbit around the 
barycenter with the relative astrometry between star and 
planet yields the mass ratio, and a dynamical mass for the 
planet. In multi-planet systems with high-precision astrom- 
etry, the the planets’ mutual gravitation can be detected by 
fitting for deviations from Keplerian motion, yielding their 
mass estimates (Lacour et al.[2021). 

The 68 Pic system is the current benchmark example of 
combining planet detection techniques to measure model- 
independent masses. Relative astrometry of 6 Pic b and 
c combined with the star's astrometric and RV measure- 
ments provides dynamical masses for the planets (Snellen 
Measurements of non-Keplerian motion due to the mutual 
gravity of 8 Pic b and c are also now possible given the 
precision of recent interferometric monitoring campaigns 
(Nowak et al.|2020). 

Combining astrometry and/or RV and direct imaging 
data can yield dynamical masses for other directly-imaged 
planets, as well as for numerous brown dwarf companions 
(Grandjean er al 2019) [Brandi e a 2020] 20216). The up- 
coming Gaia data releases for the full and extended mission 
will be used extensively in future detection and character- 
ization studies of directly-imaged planets, yielding further 
dynamical masses for directly imaged companions. The ex- 
tended mission will provide precision astrometric measure- 
ments of the host star over «10 years. 

Combining direct imaging with other detection tech- 
niques can improve planet occurrence rate measurements 
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and can better determine the planet frequency distribu- 
tionover a range of masses, orbital periods, and host star 
properties. Long-term Doppler surveys are largely com- 
plete to giant planets within ~5 au, and partially complete 
to giant planets out to ~10 au, the current inner limit of sen- 
sitivity of direct imaging surveys. Thus, combining results 
from both techniques reveals giant planet demographics out 
to ~100 au (Section|6.3). RV studies of young stars (e.g. 
allow occur 
rence rates to be compared for stars of similar ages, trac- 
ing the extent of giant planet migration over system life- 
times and its effect on planet frequency. The final Gaia 
planet catalog will also help place 
demographic measurements of young, wide-separation gi- 
ant planets into context, as astrometry will be better able to 
probe intermediate separation giant planets (~1-10 AU, = 1 
Mj,5) around younger, higher-mass stars compared to RV. 
These target stars more closely match the hosts of imaged 
planets. Combining Gaia observations of younger, higher- 
mass stars with direct imaging will illuminate trends in oc- 
currence rate as a function of stellar mass and age. 

Direct imaging also shares an important overlap with 
planet detection by microlensing. Both techniques are sen- 
sitive to planets beyond the snow line, with microlensing 
probing lower planet-star mass ratios (e.g. 
(2016), while the host mass of imaged planets can be di- 
rectly determined, and the planet mass inferred from mod- 
els. In addition, singly lensed short-period microlensing 
events can represent either a free-floating planet or a wide- 
separation bound planet (Sumi et al.[2011). Constraints on 
the wide-separation giant planet population can be deter- 
mined from imaging, thus more definitively constraining 
the free-floating planet occurrence rate. 


4. Atmospheric Characterization of Directly Imaged 
Exoplanets 


Direct imaging enables characterization of young (age 
<200 Myr) jovian exoplanets at wider separations with neg- 
ligible irradiation compared to most transiting planets. In 
this section, we first summarize from an empirical stand- 
point what we have learned about the atmospheres of the 
current cohort of young directly imaged giant planets from 
their photometry and spectroscopy. We then consider the 
theoretical side of the picture, in particular, the state-of-the- 
art in how we model these complex, cool atmospheres. 


4.1. Empirical Constraints from Time-Averaged Ob- 
servations 
4.1.1. Photometry 


Photometry in the major near-to-mid IR passbands — J, 
H, K, Ly, M, — covers the bulk of emission for young su- 
perjovian planets and provided the first empirical diagnos- 
tic of young, directly-imaged planet atmospheric properties. 
Figure |8| displays a typical near-IR color-magnitude dia- 
gram for selected directly-imaged planets and other planet- 
mass companions compared to older, more massive field 
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Fig. 8.— Near-infrared color-magnitude diagram of field M, 
L and T-dwarfs (gray circles), isolated objects exhibiting evi- 
dence of youth or possessing low surface gravity (blue squares), 
and directly imaged substellar/planetary-mass companions to stars 
and brown dwarfs (red diamonds). We annotate selected imaged 
exoplanets/planet-mass companions. Young, low-gravity objects 
form a sequence displaced redward from that of older, higher- 
gravity field objects. Data obtained from the UltracoolSheet. 


brown dwarfs and younger brown dwarfs whose ages and 
masses partially overlap with those of most imaged planets. 
Imaged planets and planet-mass companions span the full 
luminosity range characteristic of mature late-M, L, and T 
field brown dwarfs (vertical axis). However, young brown 
dwarfs, planets, and planet-mass companions typically have 
redder colors than field objects. 

The most pronounced differences between planet and 
field brown dwarf photometry occur near the transition from 
methane-poor L-type dwarfs to methane-absorbing T-type 
dwarfs. The first imaged planet-mass companion — 2M 
1207 B — and some of the first imaged exoplanets (e.g. HR 
8799 b) are particularly discrepant, appearing to populate a 
previously empty part of these diagrams consistent with a 
reddened extension of the L dwarf sequence to lower lumi- 


nosities (Chauvin er al 04] 2008 
[et al.|[2011| [Barman et al.|2011a). Over the full 1-5 um 


spectral range, photometry for these objects appears red- 
der and more blackbody-like (Figure top-left). More 
recent studies show that a number of other young imaged 
planets/planet-mass companions — e.g. HD 95086 b, TYC 
8998-760-1 c, HD 206893 b, 2M 2236+4751 B — also popu- 


late this region (e.g. 2013bj Bohn et al.|2021 
Delorme et al.[2017a][Bowler er al- POTN). 


Differences between field and young objects with T 
spectral types are less clear. Few young, T-type planet-mass 
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companions have been identified in the last decade: e.g. 51 
Eri b (Macintosh et a 2013), GU Pse B (Naud ei al 2014). 
While GU Psc B follows the sequence of field brown dwarfs 
in near-IR color-magnitude diagrams, the exoplanet 51 Eri b 
is redder than field brown dwarfs. Four additional compan- 
ions orbiting primary stars with intermediate ages (<1 Gyr; 
HN Peg B, ROSS458C, BD204-39B, Coconuts-2b, 
Xi 
show similar but less pronounced deviations with respect to 
field dwarfs in the same luminosity range. Some directly- 
imaged L/T transition planets and planet-mass companions 
have color-magnitude diagram positions discrepant from 
field objects in thermal IR passbands probing methane ab- 
sorption (3.344m, Ms) but have positions similar to field ob- 


jects in other passbands (e.g. Lp, 4.05 um) (Skemer et al. 
2012)/2014a}|Currie et al.12014a). 


4.1.2. Spectroscopy 
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Fig. 9.— Near-infrared spectra of imaged exoplanets obtained 
using low-resolution integral-field spectroscopy: from top-to- 
bottom, 8 Pic b (Chilcote et al. 2017), « And b from Sub- 
aru/CHARIS taken in a lower-resolution single-shot mode 
rie et al.|2018), HD 206893 b (Ward-Duong et al.|2021), PDS 70b 


from VLT/SPHERE (Mesa et al./2019), and 51 Eri b 
2017) from Gemini/GPI. The spectra are normalized by the flux 


between 2.05—2.15um (1.6um for 51 Eri b for clarity). Dashed 
lines indicate the offset applied to each spectrum. 


Over the past decade, IFS instruments, especially those 
used in combination with extreme AO systems (P1640, 
SPHERE, GPI, SCExAO-CHARIS) have provided critical 
low-resolution (R~20 to 80) near-IR (1-2.5,4m) spectra of 
most directly-imaged exoplanets, a representative sample of 
which is displayed in Figure[9] Near-IR spectra provide key 
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diagnostics of brown dwarf and planet atmospheres. Spec- 
tral indices derived from low-resolution data like H50 yield 
coarse estimates of spectral types (Allers and Liu|2013). 
The H-band continuum index and (possibly) the H2-K in- 
dex is a diagnostic of surface gravity 
(Canty er at 2013). 


Following trends from photometry, the spectra of young 
directly-imaged planets and other planet-mass companions 
show significant differences with the spectra of field brown 
dwarfs. Among the most notable are the following: 


e Chemistry — The L/T transition traces the onset of 
methane absorption. However, spectra of some L/T 
transition exoplanets and planet-mass companions — 
e.g. the HR 8799 planets, 2M 1207 B - show a 
lack of methane absorption compared to field dwarfs 
with similar temperatures and CMD positions (Bar-| 
man ei al. O1 Ta]b 
et al.|2018). Thermal IR spectra confirm that the HR 
8799 planets also have weaker absorption in the 3.3 


um methane filter (Doelman et al.|2022). 


Gravity — The H-band spectra for L dwarf and L/T 
transition directly imaged exoplanets / planetary- 
mass objects such as HR 8799 b, k And b, ROXs 
42Bb show highly peaked H-band spectra and/or 


red K-band spectra compared to field dwarfs (Lucas 
t al.|2001||Barman et al.|2011aj |Currie et al.|2018 
014b; |Bowler et al.|2014| [Allers and Liu|2013 


A Gauza et al.|2015). The H and K-band 
shapes probe collisionally-induced absorption (CIA) 


of hydrogen; lower gravities result in weaker CIA 
and thus peaked H-band peaks and redder K-band 
slopes. 
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Dust — At least some of the HR 8799 planets have 
spectral properties and molecular absorptions that are 
well matched by those of young free-floating objects 
at the L-T transition (see also Section (4.3). How- 
ever, no object yet reproduces the available 1-2.5um 
spectra of HR8799b, which may be due to the lack 
of identified young free-floating T-type objects. A 
handful of early-T-type objects such as the AB Do- 
radus member 2MASS J13243553+6358281 
(see Section |4.1.1) can reproduce the 
spectral bands of HR 8799 b provided that an extra 
layer of extinction by sub-micron dust particles is ap- 
plied to these empirical template spectra to match the 
spectral slope of the planet (Bonnefoy et al.|2016). 
HD 206893 b presents an even more extreme exam- 
ple, with an even flatter, more blackbody-like spec- 
trum yielding an extremely red spectrum from 1 to 
2.5 um, likely due to substantial atmospheric dust 


(Milli et al.|2017||Delorme et al.\2017b\|Ward-Duong 
et al.)2021). 


Spectra of the coolest and lowest mass imaged exoplanet 
known to date (51 Eri b) display a clear and so-far unique 
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detection of a methane absorption at 1.64m in the spectrum 
of an imaged exoplanet 
[2017], coincident with an enhancement in the K-band flux. 
Similar K-band flux enhancements have already been noted 
in young mid- to late-T dwarfs but the enhancement is par- 
ticularly extreme in the case of 51 Eri b and is due to the 
reduced collision induced absorption of Hə 
in 51 Eri b's lower-pressure atmosphere. 
Protoplanets display far more featureless, blackbody- 
like spectra. PDS 70 bc's spectra reveal an extremely red 
spectral continuum devoid of the strong water-band feature 
(1.3-1.4,m) expected given the observed luminosity of the 
planets [Mesa et aL]DOIS). Higher 
resolution K-band spectra show a lack of molecular absorp- 
tion in PDS 70 b's spectrum (Cugno et al.|2021). The cir- 
cumplanetary disks and/or cocoon surrounding each planets 
(see Fig 7 and [Christiaens er al 2019)|Benisty e al [021 
produce significant foreground IA 
et al.|2020) and may produce spectroscopic properties simi- 
lar to whatis seen in the near-infrared spectra of enshrouded 
class I protostars (e.g. |Connelley and Greene|2010). The 
near-IR spectrum of AB Aur b is reproduced by a 2000- 
2500 K blackbody but likewise lacks clear evidence for 
molecular absorption common in fully-formed substellar 
objects with similar temperatures (Currie et al.[2022). 
Aside from a few isolated cases (e.g. HR 8799 bcde), 
the spectroscopic properties of imaged exoplanets are more 
poorly constrained at wavelengths longward of 2.5um. The 
3-5 um range is particularly interesting since it contains 
both methane and !?CO absorption features diagnostic 
of carbon chemistry (CO-CH,) and cloud structures (e.g. 
[2020). The L,-band spectra of hot- 
ter young M7 to L3 companions presented in [Stone et al. 
show no significant difference with those of field 
dwarf counterparts. The L, spectrum of « And b is well 
matched with that of young, nearby free-floating brown 


dwarfs (Stone et al.|2020). A few free-floating planet ana- 
logues have high-quality 3-14.5 wm IRTE, AKARI, and 


Spitzer spectra 
mura\2012). Analysis of these data would benefit from an 
improved knowledge of the age and nature of these objects. 

While most spectra for directly imaged planets to date 
has been at resolution R <100, higher-resolution spectra 
are yielding promising results. Near-IR IFS instruments 
at medium resolving powers (R=2000-6000) fed by stan- 
dard adaptive-optics modules are available on multiple 8-m 
class telescopes (e.g. Gemini/NIFS, VLT/MUSE and SIN- 
FONI, Keck/OSIRIS). Medium-resolution IFS data enables 
innovative strategies for evaluating and subtracting the stel- 


lar halo with the spatial diversity (e.g., 2007 
Seifahrt et al.|2007; |Wilcomb et al.|2020). Most current 


medium resolution spectroscopic results are for young late- 
M and early-L young companions straddling the deuterium- 
burning boundary at p 70.5" (e.g., McElwain et al. 
et al.|2020). However, ADI observations on these instru- 
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ments (Lavigne ef al [2009] [Barman et al.[201 1a} [Meshkar 


et al.|2015) and improved data reduction strategies have 
recently allowed the extraction of high quality medium- 
resolution spectra of emblematic exoplanets HR8799 bc 


and HIP65426 b (Konopacky et al.||2013 
2015} 2021) deeply buried in the speckle noise. 


The planets’ spectra show fainter and narrower molecular 
bands (FeH, CO) and atomic lines (Na I, K D), sensitive 
to the surface gravity and primordial composition of the 
objects (C/O ratio, metallicity). They also reveal emis- 
sion lines (H,/0.656um, Pag/1.282um, BrJj/2.1661m) 


that trace active accretion (e.g., 
2021a) and that can be 
used both for characterizing the accretion processes (e.g., 
and 
detecting new planets 
POT). 


Novel data processing techniques such as “molecular 
mapping”, in combination with medium-resolution spectra, 
can also partly deblend the many ro-vibrational lines con- 
tained in planet spectra, which form the broad absorption 
bands observed at lower resolution (GPI, SPHERE). These 
molecular patterns can be extricated by cross-correlating a 
template spectrum to the planet spectrum in order to reveal 


faint molecular signatures (Konopacky et al. |2013 
let al.|2021a). “Molecular mapping”, applied to each in- 


dividual spaxel of an IFS datacube can simultaneously (i) 
detect exoplanets whose continuum emission would have 
been blurred by the speckle noise and residual stray light 
introduced by the instrument optics, (11) characterize unam- 
biguously their molecular content, and (iii) measure their 
Keplerian velocity (Hoeijmakers et al. 
2021). This method, as well as other supervised and non- 
supervised data processing techniques 
applied on such IFS data will criti- 


cally boost the detection and characterization capabilities 
on the next generation of IFS operating on the 8-m class 
telescopes (GTC-FRIDA, VLT/ERIS, VLT/SPHERE+ and 
GPI 2.0) and the ELTs (e.g.,|Houllé et al.|2021). 
High-resolution echelle spectroscopy (R>20,000) will 
go further to resolve individual molecular lines and atomic 
line profiles and provide a measurement of the objects" ro- 
tational broadening (v sin(i)). High SNR spectra have 
been collected on a few young late-M and early-L com- 
panions and free- 
floating analogues, but these observations usually require 
a large amount of telescope time. Cross-correlating lower 
SNR companion spectra or spectral residuals with templates 
can access the companion’s Keplerian and rotational veloc- 
itiess, check for specific molecular signatures, and disen- 
tangle companions from the dominant speckle noise pro- 
duced by the host star. This technique is being used on a 
steadily increasing number of companions (Snellen et al.| 


; Schwarz et al. 2018 
Xuan et al.|2020| Wang et al.|2021) and is a promis- 


ing method for identifying closer-in and fainter exoplanets 
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Fig. 10.— Demonstration of key properties inferred from direct imaging observations. (top-left) Fit to HR 8799 d 1—5 um photometry 
showing that models with thick clouds better reproduce the planet data (adapted from |Currie et al.|2011). (top-right) Keck/OSIRIS 


low-resolution H-band spectrum for HR 8799 b: a sharply-peaked H band spectrum is a signpost of low surface gravity (Barman et al.| 
2011a). (bottom) Keck/OSIRIS medium-resolution HR 8799 c spectrum, revealing lines of CO and H2O (Konopacky et al.|2013). 


(Section[7). surface gravity (Marley et al.]12012). Lower gravities yield 


. : . pressure-temperature profiles more characteristic of hotter 
4.2. Inferred Theoretical Properties from Time-Averaged dwarfs and move the depth at which carbon-based chemical 


Observations reactions are quenched deep in the atmosphere, leading to 
The empirical comparisons described in the previoustwo an elevated abundance of CO 
sections, combined with atmospheric models spanning a [ley et al./2012). In the coming years, more directly-imaged 
wide parameter space in temperature, gravity, clouds, and planet detections spanning a wider range of mass, age, and 
chemistry, provide novel constraints on the atmospheric temperature will better clarify how the clouds, chemistry, 
properties of directly-imaged exoplanets. and gravity of young jovian planets evolve (see Section[/). 
Figure [I0] summarizes our current understanding of the 
key atmospheric properties of young, directly imaged gi- 
ant planets. Young exoplanets likely are dustier and have 
thicker clouds than field brown dwarfs of the same effec- The photometric and spectroscopic observations de- 
tive temperatures (top-left panel) (Currie et al.f2011 scribed in the previous two sections are time-averages — 
Bonnefoy et al.|2016). For a fixed ob- observations either capturing short, 1-2 hour snapshots or 
served effective temperature, thicker clouds translate into averaging over multiple observations covering different 
hotter temperature profiles, meaning that the temperature rotational phases. However, directly imaged exoplanets, 
is higher at a given atmospheric pressure. The 7 = 1 sur- like their higher-mass brown dwarf cousins, are dynamic, 
face is more uniform with wavelength, resulting in a flatter, rapidly rotating objects. Measuring spectral line broaden- 
more blackbody-like spectrum (magenta line) compared to ing (v sin(i)) for field brown dwarfs and monitoring their 
models with thinner clouds (black line). This effect is most photometric variability reveals rotation periods of ~3-20 


pronounced at the L/T transition — e.g. for planets like HR hours (Zapatero Osorio et al.|2006). Recently, 
8799 bcde. The redder colors of many L-type exoplanets (2021) have even found 3 brown dwarfs with rotation 


are likely explained by clouds (e.g. 2013); periods as short as ~1 hour! Rotation measurements are 
the redder color of 51 Eri b may also be due to at least scarcer for directly imaged exoplanets and their analogues, 


partial cloudiness (Rajan et al.|2017). Planets and planet- but current measurements suggest that directly imaged ex- 
mass companions at the L/T transition also show evidence ^ oplanets are also fast rotators. Spectral-line broadening 
for non-equilibrium carbon chemistry, producing a lack of | measurements of 5 Pic b indicate a rotation period of 7- 


methane absorption in K band, flat spectra at 3-4 um and 9 hours, |Snellen et al.| (2014). Rotational periods derived 
weak emission at 5 um (Galicher et al.|2011 for 27 planetary mass objects (both as companions and free- 
2013 20142). floating) from the compilation of (20202) range 


Both clouds and chemistry are strongly influenced by from 4 to 22 hours. 


4.3. Time-Resolved Atmospheric Properties: Rotation 
Measurements and Variability 
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If an object is a rapid rotator and also has asymmet- 
ric top-of-atmosphere structure (e.g. patchy thin and thick 
clouds), it will display significant variability, appearing 
considerably brighter at some rotational phases than others. 
Large variability surveys of brown dwarfs find variability is 
common (Radigan et al 2014 
Metchev et al.|2015). Variability amplitudes extend 
up to ~20% in the near-IR. Young directly imaged giant 
planets are likely to be even more variable. From a survey 
of 30 young, low-surface gravity brown dwarfs with esti- 
mated masses «25 Mj (e.g., free-floating analogues to di- 
rectly imaged exoplanets),|Vos er al.|(2019) find that 3096 of 
low surface-gravity LO-L7.5 dwarfs are variable, compared 
to 3% of field L dwarfs. Exoplanet analogues with late-L 
spectral types appear to be particularly variable, with peak- 
to-trough amplitudes >5% in a number of notable cases 


(Biller et al.[2015] [Lew et al. 2016] [Vos et a1. [2018] 


et al [2018 2020b 2020). These 


objects have spectra that are nearly identical to those of 
the archetypical HR 8799 planets (Bonnefoy et al.|2016): 
thus, their variability properties are likely similar to those 
of these bonafide exoplanet companions. Young, planetary 
mass T dwarfs are likely more variable than their old, field 
counterparts as well. Two out of the three highest ampli- 
tude T dwarf variables known, specifically SIMP 0136 and 
2M2139, have recently been confirmed as planetary-mass 
members of the Carina Near association 
Zhang et al.|202 1b). 

Variability studies, especially spectroscopic variability 
studies, are valuable probes of the top-of-atmosphere struc- 
ture (TOA) of brown dwarfs and directly imaged exoplan- 
ets. A unique map cannot be obtained from lightcurve 
monitoring. However, indicative maps can be gener- 
ated from the spectral mapping technique, which uses 
Markov-Chain Monte Carlo techniques to fit elliptical 


spots (Kostov and Apai|2013j |Karalidi et al.|2015,|2016) 


and/or planetary-scale waves (a sinusoidal bright-dark pat- 
tern confined to a single latitudinal band, 
to multi-wavelength lightcurves. Quickly changing features 
on several field brown dwarfs may indicate that the beating 
between multiple planetary-scale waves better describes the 
lightcurves observed for field brown dwarfs than elliptical, 
starspot-like spots (Apai et al.|2017] 2021). 

Different wavelengths also probe different depths in 
brown dwarf and giant exoplanet atmospheres. Mid-IR 
data probe higher-up, cooler cloudtops; the near-IR probes 
deeper into cloud decks (e.g. |Manjavacas et al.|2021). No- 
tably, TOA structure appears to vary in between this level, 
seen as phase shifts between the near-IR and mid-IR in the 
substellar object lightcurves, where successive maxima / 
minima in near-IR are offset in time relative to their mid- 
IR counterparts. Such phase shifts have been observed in 


several brown dwarfs (Yang et al.|2016) and in at least one 
free-floating planetary mass object (Biller et al.|2018). 


The studies described in the preceding paragraph have 
mostly targeted higher mass field brown dwarfs. Next gen- 
eration telescopes and instruments such as JWST will en- 
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able similar studies for free-floating planetary mass objects 
and bright exoplanet companions such as 8 Pic b and HR 
8799 bcde. Giant exoplanet companions are likely to be 
equally variable compared to highly variable free-floating 
planetary mass objects, but are much more difficult tar- 
gets for variability studies, given the technical challenges of 
obtaining high-fidelity, highly-stable spectrophotometry for 
close companions. Two studies have attempted variability 
monitoring for the HR 8799 planets using VLT-SPHERE 
[Biller et al.|2021), reaching sensitivity 
to variability with amplitudes down to 5% for periods «10 
hours, but not yielding any detections to date. 


4.4. State-of-the-Art Modeling of Imaged Exoplanets 


4.4.1. Recent developments of 1D models 


Much of our understanding of young giant exoplanet 
atmospheres comes from the constraints the observations 
described above place on atmosphere models. Atmo- 
spheric models for directly imaged exoplanets and brown 
dwarfs share a heritage with models for very low mass stars 


Chabrier et al.|2000) [Baraffe et al. 
; Allard et al. 2015) and solar 


system planets, building a self-consistent, one-dimensional 
atmospheric structure connecting a fully convective interior 
with a radiative atmosphere (“radiative-convective” model, 
see review by [Marley and Robinson|2015). However, the 
cooler temperatures and lack of core H fusion for brown 
dwarf and exoplanet atmospheres require a number of dif- 
ferent "ingredients" compared to stars, most specifically, 
a treatment of clouds and hazes in these atmospheres and 
the ability to handle disequilibrium chemistry and vertical 
mixing. 

The very red colors observed for young planet-mass ob- 
jects cannot be described by simple 1D models without 
additional cloud or convection prescriptions 
[2012]. The most common parameterization used in model- 
ing brown dwarf and giant exoplanet clouds utilizes a sedi- 
mentation efficiency parameter fseq, which tunes cloud par- 
ticles sizes necessary to balance vertical transport of con- 
densible gases (e.g. the Eddysed grid presented in Acker-] 
[man and Marley|2001). The Sonora model grid is an update 
to the original Eddysed grid models, reaching cooler effec- 
tive temperatures and a wider range of C/O values than the 
earlier grids. To date, only the cloudless model grid is avail- 


able (Marley et al./2021), but grids incorporating a similar 


[seq cloud parameterization should be available soon. The 


ATMO 2020 models (Phillips et al.|2020) also provide a key 
update to earlier grids, in particular, the COND and DUSTY 


models of (2002), incorporating a new H-He 
equation of state, updated molecular opacity, and an im- 
proved treatment of the collisionally broadened potassium 
resonance doublet. Finally, models like Exo-REM intro- 
duce a self-consistent cloud model, with cloudy and clear 
columns, a similar f, «4 method of cloud parameterization 


(Charnay et al.|2018), and grids that consider non-solar 
compositions (Nowak et al.|2020). 
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Canonical models utilize silicate dust clouds to describe 
young exoplanet atmospheres, especially at the L/T transi- 
tion. However, other work explores alternatives. For exam- 
ple, propose that a thermo-chemical 
instability producing an extra source of diabatic convention 
can alter an atmosphere’s temperature-pressure profile, pro- 
viding another path to produce the observed red colors with- 


out clouds (see also:|7remblin et al.120171|2019||2020). 


4.4.2. Inversion techniques: achievements and challenges 


The models discussed in Section |4.4.1| are full self- 
consistent radiative-convective equilibrium models: for- 
ward models where the full 1D structure of the modeled 
atmosphere is solved iteratively (similarly as with stellar 
models). These models produce grids of synthetic planet 
spectra exploring the impact of a limited set of model free- 
parameters. The comparison of these models to observed 
spectra is known as "forward modelling". 

Bayesian inference methods such as the Markov-Chain 
Monte-Carlo method or the Nested sampling algorithm can 
be coupled to precomputed grids of forward models to in- 
terpret spectra. The method can account for correlated and 
uncorrelated noise in the data and can estimate posterior 
probability distributions on the model parameters and ev- 
idence degeneracies between those parameters, thus iden- 
tifying the sub-families of models that best represent the 


data (Rice et al.[2010 Samland et al.|2017] |Stolker et al. 
2020} Zhang et al.|2021d). Recently proposed random for- 


est algorithms are a more computationally rapid alternative 
to the Bayesian framework and can provide information re- 
garding the degree in which datapoints at given wavelengths 
constrain a specific free-parameter of the model (e. g.,{Ore-| 
shenko et al.\2020). 

Modern self-consistent forward models track both the ef- 
fect of cloud condensates on the sequestration of atomic 
species in the gas phase and the consequences of non- 
equilibrium carbon chemistry, if present. If not prop- 
erly treated, the presence of clouds and non-equilibrium 
chemistry can lead to biased determination of bulk at- 
mospheric compositions. Nonetheless, for self-consistent 
models that carefully take into account these processes, 
"forward modelling" can provide accurate fits to medium- 
resolution spectra of companions consisting of several thou- 
sands of datapoints (see Figure[I 1]. It however requires re- 
interpolating the grids of models and relies on a number of 
pre-determined sets of physical ingredients in the models. 

The complementary "retrieval" approach starts from the 
spectra and works backwards. In this case, a parameter- 
ized pressure-temperature profile is adopted and other fun- 
damental parameters (mass, effective temperature, cloud 
properties, abundances, etc.) are then retrieved given the 
observed spectrum by using Bayesian or random forest 


techniques (Madhusudhan and Seager| 2009). Retrievals 


have been used extensively to interpret transiting exoplanet 


spectra (see. Madhusudhan|2018| for a recent review). In the 


last few years, retreivals have also been applied to direct 
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spectroscopy of brown dwarfs (Line et al. [2017 
2017||2021) and extrasolar giant planets (Lavie| 
2017 2020). However, retrieval ap- 


proaches lose clear physical interpretations of the spectra. 
This approach is still in its infancy and will gradually bene- 
fit from improvements made to forward models. 
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Fig. 11.— Keck/OSIRIS medium-resolution (R~4000) K-band 
spectrum of HR8799 b from|Barman et al.| (2015) fitted with the 


Bayesian forward modeling code ForMoSA (Petrus et al.|2020 
2021) and the most recent self-consistent 1D models Exo-REM 


(Charnay et al.|2018 2021), exploring different abun- 


dance values (C/O, M/H). The fit provides a consistent representa- 
tion of the spectral continuum and !?CO absorption, confirms the 


solar C/O ratio of the planet inferred in[Barman et al.| (2015) from 


the same data, and matches well the C/O value measured for the 


primary star (Wang et al.|2020). 


4.4.3. The need for 3D models 


The majority of modeling approaches for directly im- 
aged exoplanet atmospheres are one-dimensional. How- 
ever, high-amplitude, quasi-periodic variability detected in 
both brown dwarfs and free-floating planetary mass objects 
indicate significant asymmetric top-of-atmosphere struc- 
ture, often attributed to patchy thin and thick cloud cover 
(Apai et al.|2013). (See Section [4.3] for additional details 
on these observations.) Coupled with their relatively rapid 
rotation rates which likely drive significant winds, these are 
dynamic and quickly varying atmospheres — and certainly 
not one-dimensional. Full 3D general circulation modeling 
(GCM) including a treatment of radiative transport is still 
on the horizon for brown dwarfs and directly imaged gi- 
ant exoplanets, but a number of groups have implemented 
modeling approaches that take into account the significantly 
asymmetric top-of-atmosphere structure and strong atmo- 
spheric circulation of these objects. (2010) 
use the COSBOLD code to perform local, 2D radiation hy- 
drodynamics simulations, incorporating dust clouds in the 
atmospheres of late M stars and brown dwarfs and find- 
ing that convectively excited gravity waves are an essen- 
tial mixing process in these atmospheres. For very cool 
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T and Y dwarf atmospheres with patchy salt and sulfide 


clouds, Morley et al.\(2014) introduce both clear and cloudy 


columns in their models along with a cloud-covering frac- 
tion for each model atmosphere. 
(2013) present three-dimensional, global, numerical simu- 
lations of convection in the interior of brown dwarfs and di- 
rectly imaged giant planets. These authors consider how the 
temperature perturbations driven by this convection could 
drive top-of-atmosphere asymmetry, in particular, allowing 
patchy clouds to form near the photosphere of L-T tran- 
sition brown dwarfs. use a 
shallow water model to model the atmospheric circulation 
driven by the rapid rotation of brown dwarfs and directly 
imaged giant planets. Depending on the strength of the 
heat flux and radiative dissipation, they find two circulation 
regimes: 1) with strong internal heat flux and weak radiative 
dissipation, banded east-west jets will spontaneously form 
and 2) with weak internal heat flux and/or strong radiative 
dissipation, vortices will instead form, potentially produc- 
ing some of the "spotted" features found on brown dwarfs. 
To constrain the effect of latent heating due to condensa- 
tion on both cloud growth and top-of-atmosphere asymmet- 
ric structure, [Tan and Showman|(2017) applied a idealized 
general circulation model (GCM) including a condensation 
cycle for silicate vapor. They found that simulations with 
conditions appropriate for T dwarfs developed both local- 
ized storms and east-west jets and that circulation driven by 
the latent heating due to the condensation of silicate clouds 
can generate large-scale cloud patchiness. In a series of 
recent papers, Tan and Showman have investigated a sec- 
ond possible driver of short-time evolution of clouds and 
thermal structures in brown dwarf, specifically, radiative 


cloud feedback, first within a 1D model (Tan and Showman 
, then extending to a local 3D model (Tan and Show- 
man|2021a), and finally extending to a full global geometry 


model (Tan and Showman|2021b). 


4.5. Connecting atmospheric chemical abundances to 
exoplanet formation and migration histories 


As seen in Sections [4.1.2] and medium and 
high-resolution spectroscopy, coupled with state-of-the-art 
forward- and retrieval modeling approaches can measure 
atmospheric chemical abundances for exoplanets, in partic- 
ular, key tracer ratios such as C/O and !?CO / !3CO. The 
composition of a given planet-forming disk will vary as a 
function of radius, as various molecules (water, CO, CO2) 
condense out into ices and form snow lines at progressively 
greater distances in the disk. Thus, an exoplanet atmo- 
sphere’s chemical composition will carry the imprint of 
where in the circumstellar disk it first formed and through 
which parts of the disk it migrated before reaching its final 
orbital configuration (e.g. 
Madhusudhan|2019). The nature of this 


imprint is challenging to describe end-to-end, from the pro- 
tostellar cloud to the final observed exoplanet atmosphere. 
This requires a detailed understanding of the balance of 
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which species remain in gas phase vs. which condense into 
solids, as well as the assumption that the final atmospheric 
gas composition of the planet will reflect the disk gas com- 
position. 

made the first prediction for the at- 
mospheric chemical signatures for different formation his- 
tories. For the canonical core accretion model, in which 
a several Earth-mass rocky core is built up by the accu- 
mulation of dust-to-boulder sized planetesimals and then 
accretes disk gas and residual planetesimalg?] jovian plan- 
ets should have enhanced C/O~1 if they accrete most of 
their atmospheres from gas in the disk outside of the wa- 
ter snowline. Using high-resolution Keck-OSIRIS spec- 
troscopy, measured a super-stellar 
C/O ratio for HR 8799 c, consistent with the 
predictions for the core accretion scenario. Since this first 
measurement, C/O ratios have now been published for all 
four HR 8799 planets, 8 Pic b, HIP 65426 b, and & And b 


Ruffio et al.|2021). Discrepancies in values be- 
tween early measurements and more recent measurements 
(e.g. between|Lavie et al.[2017] [Wang et al.|2020) likely oc- 
cur because higher spectral resolution measurements yield 
more accurate C/O ratio measurements. 

Interpreting C/O ratio measurements requires a careful 
comparison to disk/formation model predictions relative to 
different snowlines. To fully interpret C/O ratio measure- 
ments, several of the observational studies considered above 
also incorporate significant disk modeling in their analysis 


(cf.|Molliére et al.12020; Gravity Collaboration et al.\2020 
[Petrus et al.|2021). The details of these comparisons are 
beyond the scope of this article — see |Madhusudhan|(2019) 


for a recent review. 

A number of other elemental or isotopic ratios may yield 
insight into the formation mechanisms / locations of young, 
giant exoplanets. suggest 
that the abundances of nitrogen-bearing species may also 
trace the formation location of a given planet in the disk 


(Oberg and Wordsworth|2019). The !?CO / 1? CO ratio may 


also distinguish between "brown-dwarf" like and “planet- 


like" formation scenarios (Zhang et al.|202 lalc). 


5. Architectures of Directly Imaged Planetary Systems 
5.1. 


Given their long orbital periods, directly imaged planets 
are not currently amenable to rapid orbital characterization. 
At least one years’ worth of astrometric measurements are 
typically required in order to place initial constraints on the 
parameters of the visual orbits of planets that have been de- 
tected thus far (e.g., [Pearce et al.|2015] Blunt et al.|2017). 
In fact, no planet yet discovered via direct imaging has a 
“closed” orbit, i.e. with relative astrometric measurements 
covering the full phase of the orbit. Despite the challenges 


Orbits and Dynamics 


5Planet formation models are discussed extensively in Section|6.4 
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Fig. 12.— Visual orbits of the four planets in the HR 8799 sys- 
tem consistent with the astrometric measurements and assuming 


a dynamically stable co-planar configuration (Wang et al.12018). 


These orbits are a small subset of all visual orbits, stable or other- 
wise, that are consistent with the data. 
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associated with the analyses of these incomplete datasets, 
a great deal of information can still be learned regarding 
the orbits of imaged exoplanets, and the architectures of the 
planetary systems in which they reside. 

The HR 8799 system is probably the best-studied system 
in terms of simulations of its dynamical configuration and 


long-term stability (Fabrycky and Murray-Clay 010 
Sudol and Haghighipour 2017 
Wang et aL 018). The four plan- 


ets’ astrometric measurements are consistent with a co- 
planar configuration in, or close to, a 1:2:4:8 period ra- 


tio (Figure [12] |Gozdziewski and Migaszewski|2018] 
2018). Integrating the system forward (or backward) 


demonstrates that such a configuration is dynamically sta- 
ble, at least for many multiples of the current age of the 
system (GoZdziewski and Migaszewskij2018). While a res- 
onance between the four planets is not necessary for long- 
term stability, resonances between pairs of the inner plan- 
ets are required if their masses are = 6Mjup 
[2018). In addition to testing the dynamical stability of the 
system, these simulations have also been used to investigate 
the parameter space in which additional planets could ex- 
ist in a stable configuration 
[2018). either interior to the innermost known planet or in 
a wide orbit dynamically sculpting the outer debris disk 
(Booth et al. [2016] Read er al. 2018). 

The dynamical stability of the other known multi-planet 
systems (8 Pic, PDS 70, and TYC 8998) have not yet been 
studied at the same level of detail. Their astrometric mea- 
surements are consistent with a dynamically stable config- 
uration given the masses derived from either evolutionary 


models (Wang et al.|2021; 2020b), or from com- 


bined astrometric and radial velocity measurements 
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2020 2021b). Continued moni- 


toring of these systems is warranted to investigate not only 
their dynamical stability given revisions to their orbits and 
estimated masses, but also the influence of the planets on 
the structure of any resolved circumstellar material in each 
system (see Section [5.2). 

Early astrometry of ( Pic b quickly revealed a near-edge 
on orbit for the planet configuration (Lagrange et al.|2010), 
consistent with an edge-on inclination of the resolved de- 
bris disk (Smith and Terrile|1984). While geometrical argu- 
ments show that a long-period planet like ( Pic b will likely 
not transit, a photometric event measured in 1981 showed a 


significant variation in the star's brightness 
[Etangs et al.{1997), consistent with an eclipse of the star by 
the planet or its Hill sphere based on the ephemeris avail- 
able at the time (Chauvin et al.[2012). Continued astromet- 
ric monitoring ruled out a planetary transit but demonstrated 
that the Hill sphere would transit the star in mid-to-late 2017 
(Wang et al.[2016). An extensive photometric and spectro- 
Scopic campaign ensued to measure, or place stringent lim- 
its on, circumplanetary material as the Hill sphere transited. 
No significant event was detected: the 1981 event was prob- 
ably not linked to the transit of an extended circumplanetary 
disk (Mol Lous ef al 2018 
thy et al.|2021). Despite this null detection, exo-cometary 
material not within the 8 Pic b Hill sphere has been detected 
transiting the star with TESS (Zieba et al.|2019| Lecavleier 
et al 2022 in press), suggesting that the earlier photometric 
variability could also be due to material not associated with 
the planet. 

Monitoring efforts are also ongoing for the the orbits 
of planets within single-planet systems (e.g., Maire et al. 
[2019] [Kammerer et al.|2021). In addition to providing con- 
straints on the basic orbit geometry, these monitoring efforts 
offer the potential for future dynamical mass measurements 
when combined with absolute astrometry or radial velocity 
measurements of the host star or planet (see Section B.5). 
The orbit geometries of imaged planets within single-planet 
systems have been used to infer the presence of additional 
companions. compared the predicted 
range of periastron and apastron distances of HD 95086 b 
to the radial extent of the two circumstellar debris disks in 


this system between which the planet resides 
2013; [Su et al.|2015). Their analysis suggested that HD 


95086 b does not have a sufficiently large eccentricity to be 
responsible for truncating the outer edge of the inner debris 
disk. Instead, an additional massive planet is required on a 
shorter orbital period, one that lies below the current sensi- 
tivity limits for this star. In other systems, a significant non- 
Zero eccentricity may provide evidence for interior massive 
companions due to their effect on the relative astrometry 
between the star and the outer imaged planet (e.g., 
et al./2019). Continued astrometric monitoring will differ- 
entiate between a true non-zero eccentricity and a spurious 
detection caused by the perturbation of an inner companion. 

Comparing the orbital geometry of planets to those of 
brown dwarfs and stars may can probe of the formation 
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MWC 758 HD 61005 foie AUK ME 


Fig. 13.— The protoplanetary disk around MWC 758 (left; SCExAO team, unpublished) and the debris disk around HD 61005 (right; 
NASA/ESA/G. Schneider) both have morphologies that can be explained by the gravitational influence of an orbiting planet. The spiral 
arms around MWC 758 can be generated by an external orbiting companion (e.g., c.f. , while the 
swept-back morphology of the HD 61005 debris disk can be explained by an inner eccentric planet (Esposito et al.|2016). 


mechanisms. For individual objects, (2020b) although some structure (e.g. spirals) could instead be due 
has used the measured misalignment between the planet's to gravitational instabilities. 


orbit, rotation and stellar spin axes as an empirical test of Imaged debris disks may show sharp edges and pericen- 
various formation mechanisms and dynamical processes. ter offsets caused by the dynamical perturbations from a 
Differences between planet and brown dwarf eccentricity — jovian planet located interior (Kalas et al.[2005). Modeling 
distributions may also probe formation mechanisms (see the sharp edges of disks may help constrain the masses of 
Section[6.4] [Bowler et al. [2020]. planets (e.g. {Kalas et al.[2008). Warps in debris disks are 
also likely planetary signposts: modeling the warps in some 

52. Planet-Disk Interactions cases can constrain planet masses (Lagrange et al.[2009). 
The many main sequence stars with IR excess found in Finally, mm imaging of highly structured transitional 
the 1980s to 2000s by the IRAS, ISO, and Spitzer mis- disks may reveal indirect kinematic detection of planets or 
sions along with early resolved images of disks provided — circumplanetary material. image a cir- 
evidence for numerous, nearby optically thin, Kuiper belt- | cumplanetary disk around PDS 70 c. ALMA data for HD 
like debris disks that may be nascent planetary systems 163296 and HD 97048 may show an indirect kinematic de- 
(e.g. Smith and Terrile|1984). tection of protoplanets that currently elude direct detection 
Pluto-sized objects in these exo-Kuiper belts stir collisions Via local deviations from Keplerian rotation (kinks) 


between planetesimals to produce copious debris; stirring — |et al.|2018}|Teague et al.|2018}) Pinte et al.|2019). 
may also be driven by massive jovian planets (Kenyon and. 


2010). These 6. Direct Imaging Surveys and Population Statistics 
same missions showed evidence for numerous transitional Most of the known directly imaged exoplanets were dis- 
protoplanetary disks whose spectral energy distributions covered as part of large AO surveys. New planet detections 
(SEDs) suggested cleared cavities or gaps (e.g. and characterization of existing planets via photometry and 
, See Sect. 3.2). spectroscopy is only one part of the yield of such surveys; 
Resolved images of debris and transitional disks reveal current surveys cover samples of 500-600 stars down to ex- 
tell-tale signposts of planets (Figure [13). Images of tran- cellent contrasts and yield precise information as to where 
sitional disks resolve cavities inferred from SEDs and/or giant exoplanets reside in planetary systems. While the new 
show spiral density waves plausibly due to massive jO- — detection yield of these surveys has been relatively small, 


vian companions. For example, the disk around MWC 758 — (he deep contrasts probed enable robust studies of the fre- 
shows a complex two-armed spiral structure extending from quency of wide (>20 au) young giant planets. 


~ 25 au to over 100 au (Benisty et al.|2015); AB Aur's 

disk shows significant spiral structure on both tens of au 6.1. Survey Target Selection 

200 390 in scales, Some O w emay De dayen Dy BE Direct imaging exoplanet surveys to date have focused 
ou} primarily on stars drawn from the following samples: 1) 
[et aL [2022). ZH EDS 7O: Miele hasa cavity Ko cleared Young Moving Groups (MG): as discussed in Section 3 
by the SU two protoplanets (e.g. these groups of stars share similar kinematics and igi 


pep Recent wore suggests that actively- nated in the same star-forming region. The ongoing Gaia 
forming planets could be responsible for gaps, rings, and 


spirals seen in transitional disks (e.g. 2015alb), 


mission is helping fill in the low mass stellar members of 
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Fig. 14.— Selected high-contrast imaging results for systems in 
the Sco-Cen Association: (top) The directly-imaged wide separa- 
tion companions TYC 8998-760-1 bc (~160 and 320 au; denoted 
by arrows), and (bottom) the HD 106906 debris disk 


2021} Lagrange et al.|2016). HD 106906 also os a 27 wider- 
separation planet-mass companion (Bailey et al.|2014). 


these groups (e.g.|Gagné et al.|2018b||2021). Some notable 


moving groups that have been surveyed in depth by direct 
imaging surveys include: the 6 Pic MG, the AB Dor MG, 
the Tucana-Horologium MG, and the Sco-Cen association. 

2) Young Star-forming regions (SFR): these are re- 
gions with active ongoing star-formation and probe the ear- 


liest stages of planet formation («10 M yr) (see 


2008a|b). Notable SFRs targeted for direct imaging in- 
clude Taurus, Ophiuchus, Chamaeleon, Perseus, and Ser- 


pens. SFRs are excellent regions to search for protoplan- 
ets. They have generally been less surveyed than MGs due 
to their greater distances and confusion between exoplanet 
and disk signals. There is also some overlap between SFR 
and MG designations, with boundary cases such as TW Hya 
E Eius us within m d Sco-Cen association 


Pecaut er al 012). 


3) m i5 niche: B imaging programs fo- 
cusing on the 14+ 3 Myr Scorpius- 
Centaurus association have proven particularly fruitful, 
yielding a much higher rate of exoplanet detections com- 
pared to other young moving group surveys. At a distance 
of ~100-150 pc, Sco-Cen is the closest OB association 
(massive star-forming region) to the Sun and contains over 
one thousand young stars, well over 100 of which show 
some indirect evidence for active planet formation 
[man and Mamajek[2012). It is intermediate in properties 
between 10-100 Myr-old moving groups and «10 Myr-old 
embedded young star-forming regions. 

Figure [14] shows some direct imaging results from this 
region. Sco-Cen includes multiple directly imaged planets 
and protoplanets orbiting less than ~300 au from Sco-Cen 
host stars: T PA b (Rameau et al.|2013b), "S 65426 b 

Chauvin et al. 017), PDS 70 be (Keppler er al. 2018] Haf] 
orsa Bord. and YSES 2b ( UR END et al. ss The re- 
gion a =e a large SO of other Le eid 
planet-mass companions or objects with masses near the 
deuterium burning limit and mass ratios well below q ~ 0.1 
(ce Bohn et al D020) Bailey er al 014] TYC 8998-760- 
1 be and HD 106906 b). Follow-up observations of many 
Sco-Cen members stars with ground-based extreme AO and 
ALMA spatially resolve numerous debris disks and transi- 


tional disks (e.g. |Lagrange et al.|2016 2016 
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2016). The disks often show structures (spi- 


rals, gaps, arc) that can indicate the presence of formed or 
forming planets (see Sect. 5.2). 

4) Debris disk-bearing stars: debris disks may be sign- 
posts of imageable giant planets (see Section 5). The first 
two directly-imaged planetary systems — HR 8799 and 8 
Pic — both have luminous debris disks. Debris disk-selected 
surveys like Rameau er al 20132) 
target young, nearby stars whose debris disks have a high 
fractional luminosity (e.g. 1074-1073 

5) A and B stars: these massive stars have short life- 
times, so any companion to them must by definition be 
young («200 Myr). Statistical results from several sur- 
veys show that substellar and exoplanet companions to A 
and B stars are indeed more common than similar compan- 


ions to lower mass FGKM stars (e.g. [Nielsen et al.|2019 
see also|6.3). 


6) Accelerating Stars: unlike the categories of stars dis- 
cussed above which are selected primarily due to signs of 
youth, this category of stars show indirect evidence for the 
dynamical pull from an unseen companion. While such 
targeted surveys cannot provide unbiased demographic in- 
vestigations of exoplanets, they can produce significantly 
higher yields than blind surveys (e.g. [Crepp et al.|2016). 
Some surveys, specifically TRENDS, focused on stars with 
long-term RV trends indicating the presence of an unseen 
companion (e.g. (2018), discovering a 
number of brown dwarfs, low-mass stars, and white dwarf 
companions. A key challenge with this approach is that 
stars with sufficiently low chromospheric activity to yield 
good RV precision tend to be old and thus poorer direct 
imaging targets. 

More recently, several groups have begun using astro- 
metric data to select direct imaging targets. This method 
is immune from the age selection bias that plagues RV- 
selected target lists. For an isolated star, the proper motion 
values measured by Hipparcos and Gaia should be identical. 
However, the tug of an unseen faint companion will cause 
the position of the primary to diverge slightly from what 
would be expected for an isolated star (modulated by the 
orbital period of the companion). Differences in proper mo- 
tions measured between the Hipparcos and Gaia missions 
thus reveal accelerations pasti caused by unseen plane- 
tary d [Kervella et al.[2019] [Brandt|2018| 2021] 
Fontanive et al.|2019). While in its infr this approach is 
promising, or revealing brown dwarf companions and 
a handful of white dwarfs 
[Bowler et al.[2021). It is capable of identify- 


ing hidden planetary companions: HR 8799, 5 Pic, and HD 
206893 all have proper motion anomalies caused by known 


planetary a (Brandt et al.2021b} [Nielsen et al.| 
2019; 2021). The modest astrometric pre- 


cision of Hipparcos Eoi the current inventory of stars with 


6While (Meshkat et al.|2 find evidence of a correlation between im- 


aged planets and debris disks, debris disks may not correlate with planets 


detected by RV Yelverton et al. et al. 2020] 
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detectable accelerations likely caused by planets. However, 
future Gaia data releases should significantly expand the 
list of known stars whose accelerations may be caused by 
imageable planets. 


6.2. Summary of Surveys 


Since the advent of ground-based AO-enabled imaging 
and space-based imaging with HST, there have been several 
generations of direct-imaging exoplanet surveys: 


1. Generation 1: AO-enabled and HST efforts, 2000- 
2005 (i.e. the state-of-the-art at Protostars and Plan- 
ets V, the last Protostar and Planets conference with 
a direct imaging chapter). Typical sample size: <50 
stars. Typical contrast: A(mag) ~ 8 at 0.5”. Ex- 

(VD) 

(ESO 3.6 m, Tuc/Hor moving group), 

Lowrance et al. (HST). 


2. Generation 2: 8-m class telescope + AO + speckle 
suppression techniques (but no coronagraphy). Typ- 
ical sample size: 50-100 stars. Typical contrast: 
A(mag) ~ 10 at 0.5". Examples: 
(VLT and MMT), (Gemini), 
(VLT). 


3. Generation 3: 8-m class telescope + AO + speckle 
suppression techniques + coronagraphy. Typi- 
cal sample size: 100-300 stars. Typical contrast: 
A(mag) ~ 11 — 12 at 0.5”. Examples: NaCo 
Large Programme: : 
(2015): . NICI Sci- 


fr al] 2016), SEEDS: [Tamura] 2016); [Brandt er al. 
(2014); Janson er al.) 2013), PALMS: Bowler er al| 
2015), LEECH: 2014) 
[Maire et al (2015):[Stone et al) C018) 


4. Generation 4: 8-m class telescope + extreme AO + 
sophisticated speckle suppression techniques + coro- 
nagraphy. Typical sample size: 500-600 stars. Typ- 
ical contrast: A(mag) ~ 13 — 14 at 0.5”, using 
custom-built extreme-AO coronagraphs with spec- 
troscopic capability (due to their IFS detectors). Ex- 
amples: Gemini-GPIES (Nielsen et al.| and 


SPHERE-SHINE (Vigan er al.[2021] 
2021} Langlois et al. ). Several smaller, tar- 


geted surveys are ongoing alongside the large Gen- 
eration 4 surveys — such as the Search for Planets 
Orbiting Two Stars SPOTS (Thalmann et al.|2014. 
2018). 
the B-star Exoplanet Abundance Study (BEAST) 
2021), the Young Suns Exo- 
planet Survey SES (Bohr erar DOSS 2021), 
the SCExAO/CHARIS HGCA survey (Currie et al.| 
[2020b), and the Imaging Survey for Planets around 


Young stars (ISPY) (Launhardt et al.|2020). 


23 


Direct Imaging and Spectroscopy of Extrasolar Planets 


In aggregate, these surveys typically detect ^1 sub- 
stellar or exoplanet companion per 100 stars (e.g. meta- 
analyses by [Bowler and Nielsen|2018), al- 
though targeted surveys (e.g. focused on accelerating stars) 
may produce higher yields. 

Most of the direct imaging surveys described above have 
focused on planet detection in the near-IR, especially near 
H band. However, some surveys — e.g. LEECH and ISPY 
2020) — have searched in the thermal IR, in particular the 
Ly passband (3.78 um). For a given residual wavefront 
error, the Strehl ratio is significantly higher in the thermal 
IR than in the near IR. As shown in Figure [3] contrasts are 
more advantageous in the thermal IR for both young and 
(especially) intermediate-aged jovian planets. The planets 
B Pic b and HR 8799 e were chiefly discovered from ther- 
mal IR data; discovery papers for other planets (e.g. PDS 
70 b) also included thermal IR imaging (Lagrange et al.| 
Keppler et al. [2018). The ther- 
mal IR has lower background star contamination compared 
to the near-IR. Extreme AO systems can detect planets 100 
times fainter than conventional AO counterparts in the near- 
IR (Macintosh et al.|2014). However, thus far, extreme AO 
in the near IR has a smaller (factor of ~ a few) advantage 
over a conventional AO system operating in the thermal IR 
(e.g. see GPI/H band vs. Keck/thermal IR detections of 51 


Eri b; [Macintosh er al.[015). 


6.3. Occurrence Rates and Demographics 


The planet occurrence rate (number of planets per star) — 
and how that occurrence rate varies as a function of stellar 
or planetary properties (e.g. host mass, planet mass, orbital 
semi-major axis, host star metallicity) — strongly constrains 
models of planet formation and evolution, particularly as 
we build up larger sample sizes. Early demographics stud- 
ies focused on RV survey results, showing correlations in 
giant exoplanet occurrence with stellar metallicity (Fischer! 
and Valenti|2005) and host star mass 
and constraining occurrence rates as a function of planet 
mass and semi-major axis (Cumming et al.|2008). In partic- 
ular, [Cumming et al.|(2008) found a power-law dependence 


of planet period that corresponded to a rising occurrence 
rate of giant planets as a function of log semi-major axis 
out to ~2.5 au. 

Given the low yield of exoplanets from early direct imag- 
ing surveys, detailed occurrence rate measurements proved 
difficult. Instead, Generation 1 and 2 survey analyses fo- 
cused on constraining occurrence rates with assumed mass 
and semi-major axis distributions, or setting upper limits on 
the largest semi-major axis for which the [Cumming et al.| 
power law could be consistent with direct imaging 
yields (e.g. Masciadri er al. [2005 
. These early results showed that 


the occurrence rate of giant planets (Z1 M z) in wide orbits 
(210 AU) was relatively low (51096), and that the rising 


power law of (2008) cannot extend beyond 
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tens of au and still be consistent with the lower occurrence 
rates found by direct imaging surveys. 

Starting with the first imaged planet discoveries in 2008 
(Marois et al 20080} Lagrange eral 2009) it became pos- 
sible for surveys to directly measure occurrence rates, rather 
than simply placing upper limits. A planet of a given mass 
and separation should be easier to image around a low lumi- 
nosity, low-mass star than a higher luminosity, high-mass 
star. However, hosts (e.g. HR 8799, 5 Pic) of directly- 
imaged wide-separation giant planets were typically more 
massive than the Sun (~1.5—2 Mo). 

As Generation 3 surveys reached higher contrasts than 
their predecessors, giant planets with masses down to ~2 
M Jup became detectable at separations above ~50 au (be- 
yond ~1”) for nearby moving group stars. Yields from 
these surveys remained low, however, typically 0-2 plan- 
ets per survey. Given the limited number of nearby young 
stars, the same planet was often detected by multiple sur- 
veys (e.g. 8 Pictoris b, k And b). Occurrence rate analysis 
confirmed the relatively low upper limits implied by Gen- 
eration | and 2 surveys, with <10% of nearby, young stars 
hosting a wide-separation giant planet. Analysis of occur- 
rence rate as a function of stellar mass from these Gener- 
ation 3 surveys confirmed that wide-separation giant plan- 
ets do not follow an extrapolation of the [Cumming et al.| 
power laws. Wide-separation jovian planets were 
also found to be more common around higher-mass stars 


(e.g. 2016). Wide-separation giant planet oc- 


currence rates may also be higher for stars with Kuiper belt- 
like debris disks than diskless stars (Meshkat et al.|2017). 

Generation 4 surveys have significantly higher sensitiv- 
ity to giant planets between ~10-50 au compared to pre- 
vious generations of surveys, due to improved corongraphs 
and AO. These surveys had larger survey sizes, and detected 
planets at separations and contrasts not reachable with pre- 
vious instrumentation (e.g. 51 Eridani b at ~2 M; and 
713 au). Higher yields of substellar companions meant 
more detailed demographic studies became possible, with 
the first analysis of GPIES data detecting 3 brown dwarfs 
and 6 giant planets, and the first SHINE demographic anal- 
ysis finding 8 giant planets and 7 brown dwarfs (Figure[15]. 
These surveys confirmed at —3c that wide-separation giant 
planets (2-13 M; at 10-100 au) had higher occurrence rates 
around higher-mass (71.5 Mọ) stars compared to Solar- 
type stars Nielsen et al.[2019). 

These higher planet yields also allowed constraints the 
occurrence rate of wide-separation giant planets as a func- 
tion of semi-major axis, planet mass, and stellar host mass. 
Power law fits to the occurrence rates of wide-separation 
giant planets as a function of these parameters show that 
the giant planet occurrence rate at 10-100 au rises with in- 
creasing stellar mass, decreasing planet mass, and decreas- 
ing semi-major axis (Nielsen et al.[2019). Comparing the 
power-law fits for giant planets brown dwarfs suggests that 
planets and brown dwarfs are not drawn from the same pop- 
ulation at wider separations, consistent with separate for- 


mation mechanisms (Nielsen et al.|2019). When combined 
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Fig. 15.— Tongue plots showing detections and sensitivities to 
substellar companions for the GPIES (top) and SHINE (bottom) 
surveys. Unlike previous generations of surveys, these surveys 
could detect giant planets as close as ~10 au from their parent star 


10° 


(Figure 4 of|Nielsen et al.|2019| Figure 1 of|Vigan et al.|2021). 


with results from RV surveys at smaller separations, the dis- 
tribution of giant planets as a function of semi-major axis 
also appears (Figure[16) consistent with a peak at 3-10 au, 
near the snow line, with falling occurrence rates at smaller 
and larger separations. 


RV population studies like (2019) and 
(2021) are now sensitive to giant planets be- 


yond the snow line. Both studies find a peak in the giant 
planet occurrence rate near the snow line. The frequency 
of planets with masses of 5-14 M; between 10-100 au ex- 
trapolated from the [Fulton et al.) (2021) RV study is in ex- 
cellent agreement with the same measurement (in the same 
mass and separation bins) from the Gemini-GPIES survey 
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Fig. 16.— Occurrence rate of giant planets (5-13 Myup), brown 
dwarfs (13-80 Mnp), and substellar companions (5-80 Mj) from 
multiple direct direct imaging and radial velocity surveys. Gi- 
ant planets around higher-mass stars (top left panel) and around 
all stellar types (top right panel) are consistent with a peak in 
the distribution near the snow line (73-10 au), with falling oc- 
currence rates toward larger and smaller separations. Current de- 
mographics results provide robust evidence that giant planets are 
more common around higher-mass stars, and weaker evidence that 
giant planets and brown dwarfs follow different underlying distri- 


butions. (Figure 8 of|Nielsen et al.|(2019)) 


(Nielsen et al.\2019). Taken together, results from imaging 


and RV imply a giant planet occurrence rate (across all sep- 


arations) of ~20% (Fulton et al.|2021) for Sun-like stars, 


and potentially significantly higher for giant planets around 


higher-mass stars (Nielsen et al.|2019). 


6.4. Testing Formation Mechanisms with Direct Imag- 
ing Exoplanet Demographics 


The demographics of directly-imaged planets can pro- 
vide powerful constraints on their likely formation mecha- 
nism. There are two main pathways by which giant planets 
may form: 


1. Core Accretion — a rocky core accretes first from 
planetesimals, followed by a gas envelope (cf. 


1980} Pollack et al.|1996; Alibert et al.|2004 
2013). The proto-planet's position within 


the disk will eventually determine its final mass; for 
instance, Jupiter, closer to the peak density of the 
Sun's protoplanetary disk, managed to accrete a con- 
siderably larger gas envelope than Neptune. Pebble 
accretion can enhance the protoplanetary accretion 
rate, enabling the formation of massive planets far 


from their star (for a review, see Johansen and Lam- 


brechis[2017). 
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2. Gravitational Instability within a Disk — If a disk is 
sufficiently massive to experience self-gravity, parts 
of the disk can undergo fragmentation and quickly 


form giant planets (cf. 
let al.|2015). note that in most 
cases, these fragment-formed companions will grow 
further into substellar or low mass stellar compan- 


ions — only the lowest mass members of this cohort 
remain planetary. 


For very low mass primaries (e.g. brown dwarfs), 
binary-like formation can also produce planetary mass com- 
panions (Offner et al.[2010). This is the probable formation 
mechanism for the brown dwarf - planetary mass binary 
2M1207 B (Chauvin et al.|2004), for instance. 

Giant planets observed at wide separations experience 
significant migration or scattering either during their for- 
mation or immediately afterwards (Morbidelli et al.|2009! 
huber et al.|2021). Exoplanet demographics from large sur- 
veys provide critical tests as to which of these formation and 
migration mechanisms are dominant, as different formation 
mechanisms place planets in very different parts of their 
stellar system. Gravitational instability predicts giant plan- 
ets generally at wider separations, with more higher-mass 
planets than lower-mass ones, and no strong dependence on 
stellar mass. Core accretion, on the other hand, predicts 
more giant planets closer to the star and more lower-mass 
planets than higher-mass ones. Figure[17|shows a compari- 
son of simulated planet populations formed via core accre- 
tion (/da et al.|2013) against a population formed via disk 
instability (Forgan et al.|2015). Both of these studies simu- 
late multiple planetary embryos per system; after an initial 
formation phase, both model populations undergo an evo- 
lution phase where they can scatter or be scattered by other 
objects in their system. Obviously this is comparing just 
two simulation codes and a wide variety of similar popula- 
tions exist in the literature] However, while the exact pop- 
ulation simulated by each team and model considered will 
vary, the hallmarks of each formation mechanism remain. 

In the last decade, a number of large direct imaging sur- 
veys have used results from these types of population syn- 
thesis codes to inform their statistical analyses. Informed 


by theoretical results, (2011) 2012); [Rameau] 
(20 13a) define regions in mass-semi-major axis space 


where gravitational instability formation would be possible. 
Given non-detections in these mass - semi-major axis re- 
gions in their surveys, they were then able to set upper limits 
on the frequency of giant planets formed via gravitational 
instability. (2017) found that the disk instabil- 
ity population could describe the companions detected in 
their survey, but only if formation by disk instability was 
relatively rare. Comparisons of core accretion population 


"This topic is extensively covered in Chapter 20 from this volume, "Planet 


Formation Theory in the Era of ALMA and Kepler: From Pebbles to Exo- 
planets". 
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synthesis models to direct imaging surveys before the ad- 
vent of Gemini-GPI and SPHERE-SHINE (e.g. 
Vigan er al 017) did not 
yield meaningful constraints as most of the predicted core- 
accretion planets were undetectable for these surveys. 
Statistical results from the Gemini- GPIES and SPHERE- 
SHINE surveys sug- 


gest that core accretion is the dominant mechanism form- 
ing directly imaged giant planets. 
compared the power-law fit to 10-100 AU giant planets 
detected in the GPIES survey to predictions of both gravita- 
tional instability and core accretion. While the brown dwarf 
companion population was well-described by gravitational 
instability, found giant planets closer 
to their star in comparison to the brown dwarf companions, 
with more lower-mass planets than higher-mass planets, 
and with a higher frequency of planets around higher mass 
stars in the samples — all hallmarks of core accretion. Fitting 
both disk instability and core accretion populations(Forgan| 
simultaneously to the de- 
tected companions and contrast limits of a sample of FGK 
stars observed as part of the SPHERE-SHINE F150 sample, 
qualitatively found that the contribution 
of the core accretion simulated population is larger than 
the contribution of the gravitational instability simulated 
population, suggesting that core accretion is the dominant 
formation mechanism for this sample as well. 

Directly imaged giant exoplanets are numerous enough 
to study the demographics of companions themselves, 
rather than only from detection statistics across a survey. 
For instance, considered eccentrici- 
ties of 27 long-period giant planets and brown dwarf com- 
panions between 5 and 100 au, derived from their own 
orbital fits to literature and additional Keck astrometry. 
They tentatively found (~2-c) that their popuation dom- 
inated by giant planets (2-15 Mz) had somewhat lower 
orbital eccentricities than a brown dwarf-dominated pop- 
ulation (15-75 Mj), which displayed a broad distribution 
of eccentricities with a peak in eccentricity between 0.6 
and 0.9. These population-wide differences in eccentricity 
may suggest different formaiton mechanisms for the two 
populations: formation within a disk for directly-imaged 
exoplanets vs. cloud fragmentation (binary-like formation) 
for brown dwarf companions. 


7. The Future of Direct Imaging 


Over the next 5 years, new technology applied to exo- 
planet imaging instruments from the ground will close the 
gap between the planets we can image and those we have 
detected indirectly, while JWST will enable the detection 
of planets undetectable from the ground. By 2030, new fa- 
cilities should provide the first reflected-light detections and 
spectral characterizations of mature planets. Finally, by the 
2040s, we should be able to directly detect rocky, habitable 
zone planets around the nearest M stars and Sun-like stars 
using ELTs and NASA/ESA flagship missions. 
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Fig. 17.— Plot reproduced from (2017) com- 


paring simulated population of companions formed via core ac- 
cretion (blue/green 2013) and gravitational instability 
(red/yellow Fe OTE Lohr colors represents higher 
occurrence rates. Histograms of the distributions for each model 
population in semi-major axis (au) and mass (Mz) are shown on 
the top and side respectively. Mean detection probabilities and 
the detected companions from the |Vigan et al.) (2017) sample are 
overplotted as black contours and points respectively. Both of 
these simulated populations include multiple planetary embryos 
per system; both model populations undergo an evolution phase 
where they can scatter or be scattered by other objects in their sys- 
tem. Nonetheless, disk instability predominantly forms massive 
companions (1-100 M; far from the star, whereas core accretion 
dominates within 10 au and predominantly produces lower mass 
companions. 
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In this section, we project the near, medium, and long- 
term development of direct imaging through technological 
innovation. We describe as well new facilities applying these 
innovations and the science questions we will begin to an- 
swer via these new capabilities. 


7.1. 


The hardware enabling current high contrast imaging ca- 
pabilities described in 82.2.1|- detectors, deformable mir- 
rors, high performance computing — will continue to ad- 
vance and provide deeper contrast levels. Improved coro- 
nagraph designs will also open access to planets at separa- 
tions near the telescope diffraction limit. For ground-based 
extreme AO these incremental gains are being implemented 
in instruments on current large telescopes, either as part 
of a continuous development process (SCExAO, MagAO- 


X, KPIC[Guyon er al. [2020] 
2021). or major up- 
grades to existing facility instruments (GPI-2.0, SPHERE+ 
Boccaletti et al.[2020). The same hard- 
ware technologies are at the core of future ground and space 
instrument designs, and will provide science capability im- 
provements beyond the already significant gains expected 
from re-deploying current instruments on larger telescopes. 


New Tools 
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Advances in low-noise detectors and optical fiber tech- 
nologies will enable exoplanet spectroscopy at high spectral 
resolution. Early on-sky demonstration (KPIC, REACH, 
are paving the way for 
more capable future systems (MODHIS,|Dumas|2020). Po- 


larization differential imaging systems optimized for high 
contrast imaging (ZIMPOL, VAMPIRES, 
[Norris et al.|2020), when deployed on largest tele- 
scopes, will be able to reveal partially polarized starlight 
reflected by exoplanet atmospheres. 

Upcoming transformative advances in wavefront con- 
trol, calibration, and starlight suppression will yield order- 
of-magnitude level improvement in contrast, or enable new 
measurement capabilities. These include: 


* Sensor fusion and predictive control will optimize 
how real-time wavefront sensing telemetry is used to 
estimate wavefront errors. In predictive control, cur- 
rent and past measurements are combined to extrapo- 
late future optimal AO commands, reducing time lag 
errors and improving sensitivity 
(2018). By combining information from multiple sen- 
sors operating at separate wavelengths or in different 
optical configurations, the precision and accuracy of 
Raus UM should also be further improved 


Guyon et al [2020). 


PSF calibration from real-time telemetry holds the 
potential to numerically subtract starlight from sci- 


ence ee at ~ ane noise level (Gilles et al. 
Guyon et al 2020). 


Coherent modulation can provide a real-time esti- 
mate of starlight, which can then be numerically re- 
moved from science images. The modulation can be 
performed c in a self-coherent camera 


(doz et al.[2006) , or temporally with coherent dir 


ential imaging. 


Photonic technologies provide new ways to filter, 


transport and combine light (Gatkine et al.| 2019). 


Applications include fiber-fed spectroscopy, wave- 
asc sai and et interferometry ( 
Laser beacons can improve wavefront control per- 
formance beyond the fundamental limits imposed by 
starlight photon noise. The approach has been suc- 
cessfully used for non-extreme AO ground-based ob- 
servations, and is now being considered for both 


ground-based and space based (Douglas et al.12019) 


high contrast imaging applications. 


Starshades can provide high contrast imaging with- 
out requiring exquisite optical stability. They are 
particularly well-suited for long exposure wide spec- 
tral coverage spectroscopy of exoplanets with space- 


based telescopes (Gaudi et al.|2020). 
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Fig. 18.— Detectability of mature giant planets in optical re- 
flected light using an extreme AO planet-finding camera at an ELT. 
Top: simulated Gaia detected planets (drawn from the distribu- 


tion in [Perryman et al.|2014), Bottom: known RV planets. The 


contrast curve shown is a prediction for the ELT PCS planet-finder 
(Kasper et al.[2010), but is representative of the performance ex- 
pected for similar instruments on ELTs. Known RV planet proper- 
ties draw from the NASA Exoplanet Archive. Contrasts are com- 
puted at the planets' brightest phases. 


7.2. New Instruments and Facilities 


James Webb Space Telescope (2022) — Launched in De- 
cember 2021, the cryogenic, IR-optimized JWST may sig- 
nificantly expand the diversity of imaged exoplanets and 
will break new ground with exoplanet characterization. 
The 6.5m JWST will have poorer resolution and perfor- 
mance at the smallest angular separations than most leading 
exoplanet-imaging platforms. However, JWST will have 
unrivaled sensitivity in the thermal/mid IR (3-20 wm) com- 
pared to ground-based instruments, which must contend 
with bright sky emission: JWST can discover cold sub- 
Jupiter mass planets at wide separations around the nearest 
star (Carter ei aI2021). Thermal 
IR IFS data will be free of telluric contamination. 

Next-generation extreme AO on 6-10m telescopes (2022- 
2026) — Over the next few years, upgraded systems should 
improve planet-to-star contrasts by a factor of 10-100x at 
small separations, reaching contrasts floors nearing 1077, 
and improving performance for optically fainter stars. 
SCExAO on Subaru will undergo major hardware and soft- 
ware upgrades, improving upon its current implementation 
of advances such as predictive control incorporating ones 
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like real-time telemetry calibration, coherent modulation, 


and photonics technology (Guyon et al.{2020} 
2020a). MagAO-X (Males et al.|2020) will leverage recent 


hardware and algorithm advances to push high contrast per- 
formance in the optical. Upgrades to GPI and SPHERE 
will likewise feature improved wavefront sensing capabili- 
ties (e.g. [Chileote er al [020] Boccalett er a 030). 
Roman-CGI (2026) — The Coronagraphic Instrument on 
the Roman Space Telescope will provide the first test of 
focal-plane wavefront control and advanced coronagraphy 
in a space-borne environment (Spergel et al 013). CGI's 
benchmark performance is —-10 ^ ': comparable to the best 
contrasts expected from ground-based extreme AO on 8- 
10m telescopes in the next 5 years. However, laboratory 
tests simulating a space environment suggest that CGI could 
reach contrasts as deep as ~107°: sufficient to image true 


Jupiter analogues in reflected light (Shi et al.|2017] |Spergel 
[et _al.|[2013). The potential success of CGI's technolog- 


ical demonstration may justify a follow-on science pro- 
gram focused on detection and spectral characterization of 
reflected-light gas giant planets. 

Planet imaging with Extremely Large Telescopes (2030s) 
— The large apertures and excellent sky conditions at the 
sites of the upcoming European Extremely Large Telescope 
(E-ELT; 39m; Cerro Armazones), Thirty Meter Telescope 
(TMT; 30m; preferred site is Maunakea), and the Giant 
Magellan Telescope (GMT; 24.5m; Las Campanas) make 
exoplanet imaging with ELTs a critical science case. All 
three telescopes have approved or planned high-contrast 
imaging instruments covering the red optical to the mid-IR. 


Mid-IR instruments — e.g. TMT/MICHI (Packham et al. 
2018), E-ELT/METIS (Brandl et al.|2021) — use a single 


high-order AO correction plus a coronagraph to achieve ex- 
pected contrasts of ~10~ at 3-10 um capable of imaging 
warm Neptune or larger-sized planets around many nearby 
stars and exo-Earths around a dozen or so solar-type stars 


(Lopez-Morales et al.|2019). Those operating in the red op- 
tical to near-IR — e.g. TMT/PSI (Jensen-Clem et al.|2021) 
and E-ELT/PCS 
ond AO correction in series to achieve expected contrasts 
of «1075, sufficient to image rocky, habitable zone planets 
around nearby M dwarfs. 

Space-borne high-contrast imaging flagship missions 
(2040s) — NASA is considering two potential exo-Earth 
imaging instruments with apertures ranging between 4m 


and 16m: HabEx and LUVOIR (Gaudi et al.|2020 
LUVOIR Team|2019), and the recent NASA ASTRO2020 


decadal review rated large-aperture, space-based planet- 
imaging telescopes such as these as its highest priority for a 
flagship mission. Both designs will employ advanced wave- 
front control and coronagraphy, building upon the experi- 
ence of technical developments from Roman-CGI and the 
ground and yielding 101? contrasts. Some designs would 
fly with a starshade capable of delivering even deeper con- 
trasts and higher SNR detections of exo-Earths identified 
from coronagraphy. The European Space Agency is con- 
sidering a mission — LIFE - to detect and spectrally char- 
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acterize exo-Earths through high-contrast interferometry 


(Quanz et al.|2021). 


7.3. Future Science Focuses 


The next 5 years: young giant exoplanet atmospheres 
at A > 5 um with JWST — JWST will enable the first ex- 
ploration of directly imaged exoplanet atmospheres at > 
5 um. Especially past 10 jum, directly imaged exoplanet 
as well as brown dwarf atmospheres are almost completely 
unexplored; the only >10 um spectra that exist for any sub- 
stellar object are for a handful of L and T field brown dwarfs 
from the Spitzer cryogenic mission (Cushing et al.[2006). 

The ~10 um spectral region is particularly diagnos- 
tic of atmospheric properties. Silicate emission at 10 um 
may directly confirm silicate clouds in these atmospheres 
(Luna and Morley|2021), the existence of which has previ- 
ously been inferred by modeling, but not directly detected. 
Longer wavelengths may also better distinguish between 


different theoretical models, see e.g. (2012) 
and (2014a). JWST will also potentially de- 


tect thermal emission from sub-Jupiter mass planets >30 
au from their host stars and sub-Saturn mass planets >50 


au from their host stars (Beichman et al.|2010} 
2021). 


The next 5 years: ground-based advances — Next- 
generation extreme AO systems will achieve contrasts 
capable of detecting young Jupiters and intermediate- 
aged superjovian-mass planets at Jupiter-like separations. 
Should contrasts at or slightly below 107^ be reached, 
these new instruments might be capable of providing the 
first images of planets in reflected light. 

Direct imaging detections of significant numbers of 
super-Earth or Earth mass planets will have to wait for 
the advent of ELTs and new space telescopes. However, 
imaging the known habitable zone planet Proxima Cen- 


tauri b (Anglada-Escudé et al.|2016) might be possible 


with 8-m class telescopes and substantial dedicated ob- 
serving time. For example, (2017) propose 
targeting this planet by combining VLT/SPHERE with the 
ESPRESSO high resolution spectrograph. They estimate 
that 20-40 nights of observations would yield a 5-ø detec- 
tion of the planet and that 60 nights of observations would 
yield a 3.6-c detection of the biosignature gas O2. 

The next 10 years: Imaging and characterizing mature 
giant planets — With the advent of ELTs on the ground and 
dedicated space-based missions such as CGI, direct imag- 
ing of exoplanets will expand from detecting young plan- 
ets via thermal emission in the infrared to detecting ma- 
ture planets via reflected light in the optical. Indirect meth- 
ods like astrometry and RV will provide a sample of solar- 
system age giant planets which can then be imaged via their 
reflected optical light using exoplanet imaging cameras on 


ELTS such as PCS on the E-ELT (Kasper et al.|2021) and 
PSI on the TMT (Guyon et al.|2018). For example, Fig- 
ure[18|compares a predicted ELT planetfinder reflected light 


contrast curve (Kasper et al.|2010) with a simulated popu- 
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Fig. 19.— Detecting and characterizing exo-Earths with future facilities (from |Lopez-Morales et al.12019). (left) Simulated image 


of solar system around o Cen A imaged at 10 wm with an instrument like TMT/MICHI. In several hours, an Earth and Jupiter are 
visible: in 100 hours, 4 planets are visible. (right) Biosignatures in the optical and mid IR probed by a space coronagraphic mission like 


HabEx/LUVOIR and from the ground with ELTs, respectively (Schwieterman et al.| 


lation of Gaia-detected giant planets from 
and known RV planets. Extreme AO planet-finding 
cameras on ELTs will directly image and spectrally charac- 
terize dozens of solar-system age giant planets. 

These new facilities and instruments will also enable 
deeper characterization of the current cohort of young, di- 
rectly imaged exoplanets. In particular, these instruments 
will incorporate high resolution spectroscopic capability 
alongside extreme adaptive optics and state-of-the-art coro- 
nagraphs (Kasper et al.|[2021). This will enable high- 
precision spectroscopic and rotational studies, including 
the possibility of directly reconstructing top-of-atmosphere 
(TOA) cloud structure for a handful of young giant exo- 
planets [Crossjield[2014}. A particu- 
larly exciting future opportunity will be direct mapping of 
giant planet atmospheric structure via the Doppler imaging 
technique, which uses time-resolved, high-resolution spec- 
troscopy (R>40000) to produce a high SNR composite line 
profile, then uses maximum entropy methods to reconstruct 


TOA structure from that line profile (Vogt et al.{1987). To 


date, only one extremely bright brown dwarf, Luhman 16B, 


has been Doppler-mapped (Crossfield et al.12014) and only 


2-3 others are sufficiently bright for this technique with 8-m 


class telescopes (Crossfield|2014). Extremely large tele- 


scopes will enable Doppler imaging of dozens of brown 
dwarfs and a few directly imaged planets (Crossfield|2014). 

The next 20 years: Imaging and characterizing exo- 
Earth twins — ELTs will play an important role in imag- 
ing and characterizing habitable zone rocky exoplanets 
around a range of primary star masses. Instruments like 
E-ELT/PCS and TMT/PSI will image and spectrally charac- 
terize a few dozen habitable zone exoplanets around M stars 
in red optical to near-IR reflected light at wavelengths prob- 
ing key biomarkers (e.g. Og at 1.27 uim; |Lopez-Morales| 
let a.|[2019). Very deep mid-IR observations of nearby 


solar-like stars can detect thermal emission bracketing the 
~10um ozone feature in ~10-15 systems (Figure [19] left 
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2018). 


panel). 

Using a combination of advanced wavefront control, 
coronagraphy, and/or a starshade, space-based missions like 
HabEx, LUVOIR, and LIFE should be capable of imag- 
ing and spectrally characterizing exo-Earths around several 
tens of Sun-like stars. The NASA missions will split their 
exo-Earth finding science program into a shallow corona- 
graphic search and a deep, focused follow-up of candidate 
habitable zone planets conducted with coronagraphy (LU- 
VOIR) and/or a starshade (HabEx). Survey designs are still 
in development and strongly depend on the adopted value of 
NEarth — the frequency of Earth-like planets around Sunlike 
stars. Currently, values for Ngartn are extrapolated from the 


Kepler transiting planet mission (Dressing and Charbon- 
2015; |Kopparapu et al.|2018), although these Ngarth 


values describe only the frequency of small, rocky planets, 
not necessarily habitable planets. 

These missions will provide not just detections but spec- 
tra to identify biosignatures gases (and hence life) in the 
atmospheres of these planets (Figure right panel). No 
one biosignature gas by itself can clearly indicate the pres- 
ence of life, as many biosignature gases, such as O2, can 
have abiotic sources (Meadows et al.|2018). Thus, con- 
firming the presence of life on a planet requires detecting a 
combination of gases that do not appear together in equilib- 
rium state — disequilibria in this case is driven by the pres- 
ence of life (Meadows et al.[2018). To ensure coverage of 
a wide range of biosignature gases, these missions will de- 
liver spectra over a broad wavelength coverage: e.g. from 
near-ultraviolet to near-IR for LUVOIR. 

The study of exoplanets has yielded many surprising re- 
sults, such as the existence of hot Jupiters, the ubiquity of 
super-Earths, and the extremely red colors of young giant 
exoplanets, thus, caution is merited both in our predictions 
of what we may find with future instruments as well as in 
our interpretation of the data yielded by such instruments. 
Considering this, it will likely require detections of key 
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biosignature gases in a statistically significant sample of 
Earth-like exoplanets to be firmly convinced that we have 
truly detected life outside our own solar system. 

The detection of habitable zone planets around Sun-like 
stars through M stars through direct imaging will be nec- 
essary to provide our first steps towards truly providing 
a context for our own Earth. M star planets encounter a 
very different radiation environment than planets around 
FGK stars and are predicted to be tidally locked, which 
may pose significant obstacles to true habitability 
et al.|2019). To understand the frequency of life, it is then 
key to detect and characterize both M star habitable zone 
planets and planets more like our own — the only planet 
we know to date to truly be inhabited. However, detection 
and characterization of true exo-Earth twins (defined here 
as planets with similar mass as the Earth in the habitable 
zone of Sunlike stars) pose a challenge for transit and RV 
techniques. If viewed as an exoplanet, the Earth would 
only appear to transit once a year, thus requiring observa- 
tions over many years or even decades to reach sufficient 
S/N in transit spectra to search for biosignatures. While the 
sensitivity of current RV instruments such as ESPRESSO 
and HARPS is sufficient to detect the RV signal of an exo- 
Earth twin around a Sunlike star, that signal will generally 
be drowned out by noise from the stellar activity of the host 
star Langelier et al 2021). Thus, 
while many M star habitable-zone exoplanets will be dis- 
coved first via transit or RV, then characterized with direct 
imaging, direct imaging is more naturally suited for both 
the discovery and characterization of true Earth twins. 


This review thus concludes where it began, with the 
prospects for identifying, confirming, and characterizing a 
true Earth twin around a nearby Sun-like star via direct 
imaging. The first direct images of exoplanets were pub- 
lished less than 15 years ago. Since then, exoplanet di- 
rect imaging has advanced rapidly, constraining the orbits, 
masses, atmospheres, population demographics and forma- 
tion mechanisms of young giant planets. We forecast the 
rapid technical and scientific progress driven by the many 
and varied exciting new telescopes and instruments coming 
online in the next 20 years. 

With these new capabilities, direct imaging and spec- 
troscopy are highly likely to yield the first detections of 
true Earth twins and following that, the first detection of 
life outside the solar system — possibly the first substantive 
steps towards knowing whether or not we are alone in the 
universe. 
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